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ABSTRACT 


Software  developed  for  processing  Nell  Brown  Instrument  Sys- 
tem/NORDA  profiler  data  Is  documented  In  this  report.  The  soft¬ 
ware  Includes  programs  for  translating  the  profiler  data  from 
original  NBIS  format  to  engineering  units  In  UNIVAC/NAVOCEANO  FEB 
files  and  for  editing  and  correcting  the  data  subsequently.  This 
report  provides  complete  descriptions  of  the  programs  as  well  as 
operating  Information. 
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DXUMENTATION  OF  SOFTWARE  FOR  THE 
NEIL  BROWN  INSTRUMENT  SYSTEMS/NORDA  VELXITY/CTD  PROFILER 

I.  INTRXUCTION 

This  report  constitutes  both  documentation  and  user  manual  for  software  associ¬ 
ated  with  the  Nell  Brown  Instrument  Systems/NORDA  velocIty/CTD  profilers.  As  the 
software  Is  In  a  constant  state  of  development,  this  report  Is  presented  In  a  loose- 
leaf  format,  allowing  for  replacement  of  outdated  software,  correction  of  existing 
software  and  addition  of  new  software,  and  operation  Instructions. 

The  Initial  software  was  written  In  HPL  for  the  Hewlett-Packard  9825A  calcula¬ 
tor  by  Kim  Saunders  and  Henry  Perkins.  The  basic  data  collection  and  display  pro¬ 
grams  were  subsequently  modified  and  expanded  by  James  Vega  of  Computer  Sciences 
Corporation.  The  9825A  Is  too  small  and  slow  for  the  data  processing  envisioned, 
requiring  the  use  of  a  large  mainframe  or  super-mini  computer. 

Vega  wrote  a  set  of  three  programs  for  the  translation  of  profiler  data  to  en¬ 
gineering  units  In  UNIVAC  1108/NAV0/N0RDA  Fast  Easy  Binary  (FEB)  file  format,  the 
conversion  of  these  data  to  orthogonal,  geomagnetic  coordinates,  and  the  graphic 
display  of  these  data.  The  first  two  of  these  programs  were  subsequently  corrected 
and  modified  by  Saunders. 

Fred  Hamrick  of  Computer  Sciences  Corporation  wrote  two  programs  to  compute  the 
vertical  instrument  velocity  and  correct  the  vertical  velocity.  The  algorithms  for 
these  programs  were  developed  by  Saunders  and  Perkins,  first  testing  In  rough  form 
on  the  HP  9825A. 

After  the  GYRE  cruise  to  the  equatorial  Atlantic  In  November-December  1981,  It 
was  found  that  a  serious  over-ranging  problem  occurred  during  occasionally  strong 
down-swings  of  the  instrument.  The  algorithm  and  program  to  approximately  correct 
these  errors  were  developed  and  Implemented  by  Saunders. 

These  programs  constitute,  at  present  (April  1982),  the  existing  software  de¬ 
veloped  specifically  for  the  profiler.  NORDA  Code  331  has,  however,  a  growing  suite 
of  utility  programs  for  the  processing  of  FEB  files.  The  documentation  for  these 
programs  Is  not  included  in  this  report,  but  will  constitute  a  separate  report. 

A  DEC  VAX  11/750  super-mini  computer  Is  In  procurement  for  the  sea-going  data 
processing  of  the  profiler  data.  New  programs  to  appear  In  this  report  will  soon  be 
available  and  will  Include: 

e  real-time  profiler  data  aqulsltlon, 

t  acoustic  navigation  collection, 

e  real-time  filtering  and  correction, 

e  real-time  display. 

II.  SUMMARY  OF  DATA  COLLECTION  AND  PRXESSING 

The  Nell  Brown  Instrument  System/Naval  Ocean  Research  and  Development  Activity 
3-Axis  VelocIty/CTD  profiler  system  was  designed  for  the  study  of  upper  ocean  mixing 
and  variability.  It  Is  capable  of  measuring  three  components  of  velocity,  accelera¬ 
tion,  and  magnetic  field,  as  well  as  conductivity,  temperature  and  pressure.  This 
system  and  the  Initial  phases  of  testing  are  documented  In  Perkins  et  al.  (1980)  and 
Saunders  et  al.  (1981). 


The  profiler  system  consists  of  an  underwater  unit  and  a  deck  unit.  The  data 
are  sensed  by  the  underwater  unit,  digitized,  and  then  transmitted  to  the  surface 
via  an  audio  FSK  code  over  a  single  conductor  (sea  water  return)  cable.  The  signal 
Is  led  through  a  winch  with  slip-rings  to  the  deck  unit,  where  the  data  are 
reformatted.  There  are  three  signal  outputs  from  the  deck  unit: 

e  an  audio  output  for  backup  on  audio  tape, 

e  an  IEEE  488  parallel  bus  and, 

e  an  output  for  a  digital  tape  recorder  (either  a  Kennedy  or  Dlgl-Data). 

When  the  Initial  test  of  the  profiler  was  conducted,  the  physical  oceanography 
branch  did  not  own  a  9-track  tape  drive  that  could  be  directly  connected  to  the 
third  output.  The  branch  owned  an  Interface/buffer/recorder  system  that  was  compati¬ 
ble  with  the  IEEE  488  bus.  The  Initial  data  collection  program  was  therefore  design¬ 
ed  to  transfer  data  directly  from  the  deck  unit  to  the  Interface/buffer,  occasional¬ 
ly  breaking  Into  the  bus  to  obtain  “snapshots"  of  the  data  being  recorded.  The  pres¬ 
ent  version  of  the  data  collection  program  "DATALOGGER"  Incorporates  this  design. 
However,  It  was  recently  discovered  that  because  the  interface/buffer  Is  only  a 
singly  buffered  system,  about  15%  of  the  data  are  lost  when  the  buffer  Is  emptied 
during  recording.  The  solution  to  this  problem  Is  to  record  the  data  directly  to 
magnetic  tape  via  the  third  output. 

The  data  collection  program  Is  also  designed  to  produce  quasi-realtime  plots  of 
the  data  on  one  or  two  printer-plotters  and  to  display  numeric  data  on  a  Hazeltlne 
1420  CRT  terminal  at  9600  Baud. 

Once  the  data  have  been  recorded  on  digital  magnetic  tape,  they  are  processed 
on  the  UNIVAC  1108  at  NAVOCEANO  (this  will  be  changed  In  the  near  future  to  a  DEC 
VAX  11/750).  The  program  which  translates  the  data  from  the  NBIS  format  to  engineer¬ 
ing  units  Is  "TRANSCRIBE."  This  program  will  handle  data  written  directly  from  the 
deck  unit  to  magnetic  tape  or  tape  generated  by  the  IEEE  488-Interface  /  buffer  sys¬ 
tem. 

The  data  should  then  be  plotted  and  If  there  appears  to  be  any  evidence  of 
over-ranging,  the  "PREFIX"  program  should  be  run.  This  program  makes  a  close  approx¬ 
imation  to  the  velocities  when  over-ranging  Is  encountered. 

Once  this  Is  done,  either  "VFIX1"  or  "VFIXl-S"  may  be  applied  to  the  data. 
These  programs  attempt  to  correct  the  vertical  velocity  by  computing  and  subtracting 
the  Instrument  velocity  from  the  observed  vertical  velocity.  "VFIX1"  accomplishes 
this  by  Integrating  the  vertical  acceleration,  while  "VFIXl-S"  differentiates  the 
pressure  to  obtain  the  Instrument's  vertical  motion. 

Two  utility  programs  are  also  Included.  "UNORTHOG"  Is  used  mainly  for  testing 
and  debugging  purposes  for  looking  at  the  velocity  data  In  the  original  acoustic 
(non-orthogonal)  components.  "TSERPL0T2"  Is  a  general  plot  package  to  plot  any  of 
the  variables  In  the  FEB  files  against  either  cycle  number  or  time. 

III.  DATA  STRUCTURES 

Two  primary  structures  are  Involved  in  profiler  data  processing:  NBIS  raw  data 
format  and  profiler  FEB  files.  The  NBIS  format  Is  used  only  during  the  data  collec¬ 
tion  phase.  This  Is  a  highly  packed  format  prior  to  conversion  to  engineering  units. 
The  FEB  file  structure  Is  the  standard  file  format  used  In  the  Physical  Oceanography 
Branches  at  NORDA  and  NAVOCEANO  (Hal lock,  1981).  The  details  of  these  formats  are 
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given  below  for  reference.  The  IffilS  format  Is  given  In  Figure  1  (from  the  NBIS  pro¬ 
filer  manual,  with  permission). 

The  FEB  file  variables,  after  the  transcription  phase  are.  In  order  of  position 
In  the  data  array: 


No. 

Name 

Variable  (units) 

1 

PRESS 

pressure  (decibars) 

2 

STEMP 

slow  response  temperature  (deg  C) 

3 

COND 

conductivity  (mmho) 

4 

FTEMP 

fast  response  temperature  (deg  C) 

5 

VL0CI1 

velocity  component  1  (cm/sec) 

6 

VL0CI2 

velocity  component  2  (cm/sec) 

7 

VL0CI3 

velocity  component  3  (cm/ sec) 

8 

MAGU 

magnetic  component  1 

9 

MAGI2 

magnetic  component  2 

10 

MAGI  3 

magnetic  component  3 

11 

ACCI1 

acceleration  comp.  1 

12 

ACC  1 2 

acceleration  comp.  2 

13 

ACC  13 

acceleration  comp.  3 

14 

TIME 

reference  time  (dec.  days) 

15 

RELSEC 

relative  time  (sec) 

Further  processing  programs  occasionally  will  change,  correct,  or  replace  these 
variables.  For  Instance,  after  using  PREFIX,  variable  14  Is  replaced  with  an  esti¬ 
mate  of  the  Instrument's  vertical  velocity,  determined  from  the  pressure  time  deriv¬ 
ative.  It  Is  also  possible  to  extend  the  number  of  variables  In  the  FEB  file  struc¬ 
ture,  and  at  present  the  programs  VFIX1  and  VFIX1-S  extend  the  number  of  variables 
to  17. 

The  header  blocks  for  the  PROFILER  FEB  files  are  defined  as  follows: 


ADOC 

FDOC 


(31):  available  for  alphanumeric  Information, 
sample  Interval  time 

start  latitude  (dec.  deg.) 
start  longitude  (dec.  deg.) 
time  of  start  fix  (dec.  days) 
end  latitude 
end  longitude 
time  of  end  fix 
maximum  pressure  of  sensor 
cast  start  time  (dec.  days) 


cast  end  time 


magnetic  variation 
magnetic  dip 
ship's  speed  (kt) 
ship's  heading 

dry  bulb  temp  (deg.  C) 
wet  bulb  temp  (deg.  C) 
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(23)  - 

surface  temp  (deg.  C) 

24  - 

barometric  pressure  (mb) 

(25)  - 

wind  speed  (nm) 

(26)  - 

wind  direction  (compass) 

27  " 

significant  wave  height  (ft) 

(28  - 

(29)  - 

30  - 

(31)  - 

start  time  of  profile 

(day) 

(32)  - 

N 

M 

M  N 

(hour) 

(33)  - 

N 

M 

II  M 

(min) 

34  ’ 

II 

II 

H  M 

(sec) 

(35)  - 

36  - 

end 

M 

N  M 

(day) 

(37  - 

II 

II 

H  N 

(hour) 

(38)  - 

II 

II 

N  M 

(min) 

(39)  - 

II 

II 

H  l< 

(sec) 

(40)  - 

I DOC  (1)  -  end  of  profile  flag 
(21  -  cruise  number 

(3)  -  station  number  (Id  format) 

(4)  -  relative  segment  number 

(5)  -  series  sequence  number 

(6)  -  absolute  no.  of  1st  seg.  In 


Input  tape  no. 
year  of  cast 


no.  of  bad  cycles 
no.  of  profiles  In  cast 


series 
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Figure  1.  NBIS  Raw  Data  Format 


PARAMETER 

OISPLAY  UNITS 

LS  BIT  HEIGHT 

FORMAT 

Frama  Sync, 

240  or  019 

— 

Prassura  LS8 

sea  calibration 

see  cal Ibratlon 

AC 

Prassura  MS8 

• 

Tamparatura  LSB 

degree  celclus 

0.9  a  deg.  C 

AC 

Tamparatura  MSB 
Conductivity  LSB 

mmho 

0.001  mmho 

AC 

Conductivity  MS8  ’ 

Fast  Tamp.  LSB 

degree  celclus 

0.9  a  deg.  C 

AC 

Fast  Tamp.  MSB 

AC  Signs 

part  of  tamp. 

— 

SIGNS 

Velocity  X  LSB 

and  prassura 
16383-1  m/sec 

1/16383  m/sec 

OC 

Velocity  X  MSB 
Velocity  Y  LSB 

16383-1  m/sec 

1/16383  m/sec 

OC 

Velocity  Y  MSB 
Velocity  Z  LSB 

16383-1  m/sec 

1/16383  m/sec 

OC 

Velocity  Z  MSB 
Compass  X  LS8 

ratio  only 

OC 

Compass  X  MSB 
Compass  Y  LSB 

ratio  only 

. , 

OC 

Compass  Y  MSB 

Compass  Z  LSB 

ratio  only 

OC 

Compass  Z  MSB 
Acceleration  X  LSB 

g's  x  1000 

0.001  g 

OC 

Acceleration  x  MSB 
Acceleration  Y  LSB 

g's  x  1000 

0.001  g 

OC 

Acceleration  Y  MSB 
Acceleration  Z  LSB 

g's  x  1000 

0.001  g 

DC 

Acceleration  Z  MSB 
Spara-0  LSB** 

OC 

Spare-0  MM** 

Spare-1  LSB** 

— 

OC 

Spare-1  MSB** 

TOO-O  msec 

m  sec  x  10,  msec  x  1 

BCD 

TOO-1  sac/msec 

— 

sec  x  1,  msec  x  100 

BCD 

TOO-2  mln/sec 

— 

min  x  1,  sec  x  10 

BCD 

TOD-3  hr/ialn 

— 

hr  x  1,  min  x  10 

BCO 

TOD-4  day/hr 

— 

day  x  1,  hr  x  10 

BCD 

TOO-5  day 

— 

day  x  100,  day  x  10 

BCD 

•NOTE  1.  AC  Format 

LSB  ■  128,  64,  32,  16,  8,  4,  2,  1 

MSB  -  32768,  16384,  8192,  4096,  2048,  1024,  312,  296 

2.  SIGNS  Format 

1,  1,  1,  1,  1,  Fast  Tamp.,  Tamp.,  Prassura 
where  1  -  negative,  0  ■  posltlva 

3.  OC  Format 

LSB  -  32,  16,  8,  4,  2,  1,  0,  SIGN 

MSS  •  8192,  4096,  2048,  1024,  912,  296,  128,  64 

4.  BCD  Format 

«.  4,  2,  1,  8,  4,  2,  1 

••  All  bits  In  spara  Bytas  ara  sat  to  "1"S. 

Figure  1.  NBIS  Raw  Data  Format 
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APPENDIX  A:  AN  ESTIMATOR  FOR  VERTICAL  INSTRUMENT  VELOCITY 


Let  Wq  and  3g  represent  observed  vertical  components  of  velocity  and  accel¬ 
eration  and  let  Pg  be  the  observed  pressure.  Suppose  the  data  to  be  given  over 
the  time  Interval  tj,  tg 

Define  the  Instantaneous  vertical  velocity  of  the  Instrument  as 

*l(t)  «{*  Co  ao(t)  -  g  -  y  ]dt+  0  (A- 1 ) 

NOTE:  In  the  program,  g  was  set  equal  to  zero  and  a  correction,  AC0RR*9.99,  was 
applied  to  all  measured  accelerations. 

where  g  Is  the  local  gravitational  acceleration,  a  and  y  are  corrections  to  the 
observed  acceleration  and  $  Is  the  Instrument  velocity  at  t*tj.  The  values  of  o, 
&,  and  y  are  to  be  determined  In  an  optimal  way,  as  described  below.  Nominally, 
a  *  1  and  y  =  0. 

Further  define  the  mean  vertical  velocity  of  the  instrument  during  ti,  t2  on  the 
basis  of  the  correspondl ng  pressure  change: 

"1  “fgCPoUl)  ’  Po(t2)]/(t2  -  ti)  (A-2) 


Where  p  Is  the  mean  density  of  the  water. 

The  quantities  a.B.y  are  determined  by  the  conditions 


1*2 


and 


t  (w-f  -  wq)2  dt  »  minimum 


X  _ 

f  *(wi  -  w)  dt  *  0 

ll 


(A-3) 
(A— 4) 


That  Is,  wj  Is  required  to  resemble  wq  as  much  as  possible  and  also  correspond 
to  the  known  mean  Instrument  speed. 

Formally,  the  two  equations  above  constitute  a  linear  least  squares  problem  with  a 
side  condition.  This  may  be  reformulated  by  the  method  of  Lagrange: 


1  2C(w^  -  wo)2  +  x(w i  -  w)]  dt  *  minimum 
*1 


(A-5) 
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Where  A  is  a  Lagrange  multiplier,  the  value  of  which  Is  also  to  be  determined. 
Differentiating  this  expression  with  respect  to  o,$,v  and  A  In  turn  and  equating 
each  of  the  results  to  zero,  as  Is  required  to  minimize  the  expression,  results  In  a 
linear  system  of  4  equations  for  a,$,y,A.  These  have  the  form 


AX  *  B 

Where  A  =  (A-jj)  *  (  £t2 a-j j (t )  dt) 


B  =  (Bj)  »  (  {|2bj(t)  dt) 


(A-6) 

(A-7) 

(A-8) 

(A-9) 


NOTE:  In  the  program,  the  last  quantity  In  array  B,  / w,  was  set  equal  to  the  aver¬ 
age  velocity  determined  from  the  total  pressure  change  multiplied  by  the  to¬ 
tal  time. 


The  quantities  ajj,  bj  can  be  shown  to  be 


(A-10) 


where  fj 
f2 


(bj)  =  /  fl(*0  +  9f2) 

*0  +  9f2 
*2  (wq  +  9f2) 
w  +  gf2 


ft  2 

fl(t)  -Jt  a0(t)  dt 
fb2 

*2U)  t  dt  a  t  -  tj 

C-j 


(A-ll ) 
(A-12) 
(A-13) 


Note  that  even  though  (a^j)  Is  singular,  since  the  first  three  rows  are  multiples 
of  each  other,  (Aij)  is  not.  Hence,  It  can  be  inverted,  equations  A-6  solved  for 
o,$,y  and  the  result  used  In  A-l  to  find  an  optimal  estimate  of  the  Instrument  ve¬ 
locity. 


t 
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APPENDIX  B:  PROFILER  SOFTWARE 


PROGRAM: 

PURPOSE: 


MACHINE: 

LANGUAGE: 

AUTHOR: 

FILE  LOCATIONS: 
INPUT: 


OUTPUT: 

ADDITIONAL 

INFORMATION: 

PROGRAM: 

PURPOSE: 


MACHINE: 

LANGUAGE: 

AUTHORS: 

FILE  LOCATIONS: 


DATALOGGER 

To  log  data  and  maintain  a  real  time  display  of  certain  vari¬ 
ables  when  the  profiler  Is  collecting  data  at  sea.  Data  are 
logged  by  transferring  the  data  stream  directly  to  an  IDEAS  IEEE 
488/CIPHER  DATA  Interface/buffer.  The  variables  are  displayed 
on  a  Hazel tine  1420  CRT  terminal  and  plotted  on  two  HP  7242 
prl nter/ pi  otters. 

HEWLETT-PACKARD  9825A 

HPL 

Kim  David  Saunders,  Henry  T.  Perkins  and  James  J.  Vega 
Tape  V002,  File  4 

The  program  solicits  the  following: 

Cast  Number, 

Latitude  in  degrees  decimal. 

Longitude  in  degrees  decimal. 

These  are  used  only  for  display  on  the  Hazeltlne  terminal. 

During  normal  operation  of  the  program,  different  display  op¬ 
tions  may  be  Invoked  by  pressing  the  fO  function  key.  The  pro¬ 
gram  will  then  solicit  the  type  of  display  option  desired  and 
the  unit  to  which  that  type  of  display  Is  to  be  directed.  The 
display  options  are  listed  after  the  program  listing. 

The  output  consists  of  alphanumeric  listings  of  data  on  the 
Hazeltlne  screen  and  plots  of  data  on  either  or  both  printer 
pi  otters . 


TRANSCRIBE 

To  convert  profiler  data  from  raw  form  on  the  original  data 
tapes  (or  condensed  copies)  into  engineering  units  with  the 
velocity  expressed  In  orthogonal  Instrument  coordinates.  The 
converted  data  are  stored  in  FEB  file  format  for  ease  In  further 
processing. 

UN I VAC  1108 

FORTRAN  V 

James  J.  Vega,  corrected  and  modified  by  Kim  David  Saunders 

Absolute,  Relocatable,  and  Symbolic 
Elements  -  VEGA*L IB .TRANSCRIBE 


INPUT: 


OUTPUT: 

ADDITIONAL 

INFORMATION: 

PROGRAM: 

PURPOSE: 

MACHINE: 

LANGUAGE: 

AUTHOR: 

FILE  LOCATIONS: 

INPUT: 

OUTPUT: 


The  Input  parameters  are  solicited  by  the  program  for  use  In  an 
Interactive  environment. 

The  Input  raw  data  must  be  attached  to  logical  unit  10,  for  ex¬ 
ample  by  a  series  of  statements  such  as: 

asg.tsj  10,u9s,«tapeno.» 
move  10,<<nf11es-l» 

The  solicited  Input  parameters  are  summarized  below  for  use  In 
a  batch  environment.  All  Input  Is  In  free  format. 

Lines  1-3:  Alphanumeric  documentation  42  char/llne. 

Line  4:  Cruise  no.,  station  no.,  absolute  no.  of  1st  segment  of 
the  output  file,  Input  tape  no.,  year  of  cast  (all  Integer). 

Line  5:  Time  Interval  between  samples  In  seconds,  zero, 
starting  latitude  In  decimal  degrees,  starting  longitude,  time 
of  fix  a  start  of  station  In  decimal  days,  ending  latitude, 
ending  longitude,  ending  time,  maximum  pressure  of  profile  In 
decibars,  (all  floating  point) 

Line  6:  Magnetic  variation,  magnetic  dip,  ship  speed,  ship 
heading  (floating  point). 

Line  7:  Dry  bulb  temperature,  wet  bulb  temperature,  surface 
temperature,  wind  speed  In  knots,  wind  direction  In  degrees, 
significant  wave  height  In  feet  (floating  point). 

Line  8:  Cast  start  day,  hour,  minute,  second  zero,  cast  end 
day,  hour,  minute,  second  (Integer). 

Line  9:  Sequential  file  no.  (Integer) 

Line  10:  Station  Identification  No.  consisting  of  a  3  digit 
station  and  a  3  digit  sequential  cast  No. 

Line  11:  Message  level  0-9 

File  20  contains  the  output  FEB  file. 


CONVERSION 

To  convert  from  Instrument  orthogonal  to  geomagnetic  orthogonal 
coordinates. 

UNIVAC  1108 

FORTRAN  V 

James  J.  Vega,  corrected  and  modified  by  Kim  David  Saunders 

Absolute,  Relocatable,  and  Symbolic 
Elements  -  VEGA*L IB. CONVERSION 

File  10  -  Input  FEB  File. 

File  20  -  output  FEB  File. 


ADDITIONAL 

INFORMATION: 


PROGRAM; 

PURPOSE; 

MACHINE: 

LANGUAGE: 

AUTHOR: 

FILE  LOCATIONS: 
INPUT: 


OUTPUT: 


ADDITIONAL 

INFORMATION: 


PROGRAM: 


PREFIX 

To  read  a  FEB  file  containing  raw  profiler  data,  removing  points 
where  the  vertical  velocity  has  exceeded  the  limits  for  the 
Instrument. 

UN I VAC  1108 

FORTRAN  V 

Kim  David  Saunders  (April  1982) 

Absolute,  Relocatable,  and  Symbolic 
Elements  -  VEGA*L1B. PREFIX 

Line  1)  NUIN1 ,NSEG1,NSSEG1 
Line  2)  NOUT 

NUIN1  =  unit  number  of  Imput  FEB  file 
NSEG1  *  number  of  segments  to  be  read 
NSSEG1  =  number  of  first  segment 
NUOUT  =  unit  number  of  output  FEB  file 

The  output  FEB  file  has  the  same  structure  as  the  Input  file, 
with  the  exception  that  variable  14  now  contains  a  rough  approx¬ 
imation  to  the  vertical  Instrument  velocity  estimated  from  the 
time  derivative  of  the  pressure.  The  Initial  time  of  the  cast 
(in  Julian  days  has  been  replaced). 


The  NBIS  3VCT0  profiler  measures  the  current  relative  to  the 
Instrument  by  means  of  a  three  axis  acoustic  veloclmeter.  The 
operation  of  this  type  of  current  meter  Is  described  In  the 
NBIS  Acoustic  Current  Meter  manual.  The  point  of  Interest, 
here.  Is  that  the  apparent  velocity  along  any  axis  Is  propor¬ 
tional  to  the  phase  difference  of  the  two  acoustic  pulses, 
which  are.  In  turn,  proportional  to  the  true  component  of  the 
water  velocity  along  the  axis.  Thus,  In  principle,  the 
measured  velocity  should  be  mapped  onto  the  Interval  from 
about  -100  to  100  cm/sec  (nominal).  In  practice,  this  does 
not  occur  exactly,  as  when  the  phase  of  the  acoustic  signal  is 
near  -180  or  180  degrees,  the  gate  opening/closing  signals 
become  ambiguous.  This  results  In  random  output  velocities 
when  the  true  velocity  component  along  the  axis  Is  within  a 
small  "dead  band"  of  the  velocity  extremes.  This  program  Is 
designed  to  search  the  data  for  probable  occurrences  of  this 
over-ranging  and  to  correct  (as  much  as  possible)  by  substi¬ 
tuting  the  projection  of  the  vertical  Instrument  velocity, 
determined  by  the  time  derivative  of  the  pressure,  for  the 
components  along  the  ql  and  q3  axis. 

VF1X1 

Corrects  the  vertical  component  of  current  velocity  as  measured 
by  the  NORDA  3-component  profiler. 


PURPOSE: 


INPUT:  (Free  format) 

Line  1  -  Input  file  specifications 
TOUR  -  Unit  no.  for  Input 
NSEG  -  No.  of  segments  to  process 
NSSEG-  No.  of  starting  segment 
MSGR  -  Message  level  for  Input 
Line  2  -  Output  file  specifications 
IUNW  -  Unit  no.  for  output 
MSGW  -  Message  level  for  output 


NOTES: 


1)  The  Input  file  Is  presumed  to  be  In  geomagnetic  coordinates  as  produced  by 
program  CONVERSION. 

2)  Input  variables  are  Identified  by  name  as  follows: 

Vertical  current  speed,  uncorrected  -  VL0CG3 

Vertical  component  of  acceleration  -  ACCLG3 

Pressure  -  PRESS 

Time  -  RELSEC 

3)  The  output  file  has  the  same  structure  as  the  Input  file  except  that  the 
corrected  vertical  velocity,  named  W,  Is  Inserted  In  each  data  cycle  Immediately 
before  the  uncorrected  vertical  velocity,  and  the  instrument  velocity,  named  HI, 
Is  Inserted  In  each  data  cycle  Immediately  after  the  last  variable  (RELSEC). 


METHOD:  An  estimate  of  vertical  Instrument  velocity  w.  Is  found  from 

the  observed  acceleration  and  pressure.  Derivation  of  w.  Is 
given  In  the  Appendix.  The  corrected  velocity  w  Is  then  obtain¬ 
ed  from  the  observed  velocity  wQ  by 


PRINTED  OUTPUT:  For  each  output  segment,  the  following  quantities  are  listed: 

o  Start  and  end  times  (RELSEC) 
o  Start  and  end  pressure 
o  o,0, y, A  (see  Appendix  for  definitions) 

MACHINE:  UNIVAC  1108 


FILE  LOCATIONS: 

PROGRAM: 

PURPOSE: 

MACHINE: 


Absolute  element  C0DE331*FCHFILE1.VFIX1 
Mapping  element  C0DE331*FCHFILE1.MVFIX1 

TSERPL0T2 

To  plot  time  series  of  profiler  variables  either  versus  cycle 
number  or  relative  time. 

UNIVAC  1108 


LANGUAGE:  FORTRAN  V 
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AUTHOR: 


James  J.  Vega 

FILE  LOCATIONS:  Absolute,  Relocatable,  and  Symbolic 


INPUT:  There  are  two  Input  lines  for  each  subplot: 

Line  1:  IU,  IABSIS,  IP,  YMAX,  YSTP,  YMIN,  I0EC 
Line  2:  If  IABSIS-0:  NUMSEG,  IBEGIN,  CYCIN 
If  IABSIS*! :  TSTART,  TEND,  IPTIME 

where 


IU  ■  logical  unit  for  Input  FEB  file 
IABSIS*  0  for  cycle  number  plot 
-  1  for  time  plot 

IP  ■  position  of  variable  In  FEB  array 
YMAX  ■  expected  max.  of  variable 
YSTP  ■  lab ling  interval 
YMIN  -  expected  min.  of  variable 
IDEC  "  decimation  ratio 

NUMSEG*  number  of  segments  to  plot 
IBEGIN-  first  segment  to  plot 
CYCIN  "  number  of  cycles  per  inch 


OUTPUT: 

ADDITIONAL 

INFORMATION: 


PROGRAM: 

PURPOSE: 


TSTART-  start  time  In  seconds  (relative) 
TEND  "  stop  time  In  seconds 
IPTIME-  location  In  FEB  array  of  time 

File  25  -  output  intermediate  plot  file. 


1.  The  following  data  should  succeed  the  last  data  line  to 
ensure  proper  termination  of  the  program: 

99,0,0, 0.,0.,0.,0 

2.  This  program  Is  designed  to  use  34  Inch  ZETA  plotting 
paper.  Because  of  this,  the  maximum  number  of  variables  per 
plot  is  5. 

VF1X1-S 

To  read  a  FEB  file  containing  profiler  data  and  create  a  new 
FEB  file  irtilch  contains  a  corrected  vertical  velocity  computed 
from  the  pressure  derivative. 


MACHINE: 

AUTHOR: 

FILE  LOCATIONS: 


UNIVAC  1108 

Fred  Hamrick  (April  82) 


INPUT:  Line  1)  IUNR,IUNW,NSEG,NSSEG,NYAR,NUNV 

IUNR  ■  unit  number  for  Input  FEB  file 
IUNW  *  unit  number  for  output  FEB  file 
NSEG  •  number  of  segments  to  read 
NSSEG*  start  segment 
NVAR  *  variable  number  for  pressure 
NUMV  *  variable  number  for  vertical  velocity 

The  output  FEB  file  has  the  same  structure  as  the  Input  file 
with  two  additional  variables.  The  variable  WI  (Instrument 
velocity)  Is  written  as  the  last  variable*  and  the  corrected 
velocity  W*V-WI  (where  V*measured  vertical  velocity)  Is  written 
as  the  variable  Immediately  before  the  measured  vertical 
velocity. 

1)  Before  executing  this  program,  the  Input  FEB  file  should  be 
Interpolated  with  respect  to  time.  (This  time  difference  Is 
set  In  the  program  as  variable  DELT.)  The  Interpolation  may 
be  performed  by  executing  HTP*PROG.ZINTERP. 

2)  The  Instrument  velocity  Is  computed  as 

WI  -  RH0G1*DPDT 
Where  RH0G1  -  -  .9955 

and  DPDT  *  *  pressure  derivative 

Is  computed  as:  1=+K 

DPDT (J)  - _ 3  2  (i  *  P,+1) 

K(K+1 )  (2K+1 )*DELT  1=-K  J  1 

Where  Pj  *  pressure  values 
and  K  Is  set  to  8  In  the  program. 

(the  first  and  last  K  values  of  DPDT  are  set  to  0). 

PROGRAM:  UNORTHOG 

PURPOSE:  To  convert  from  Instrument  orthogonal  coordinates  to  Instrument 

acoustic  axes  coordinates.  (Velocity  only). 

MACHINE:  UNIVAC  1108 

LANGUAGE:  FORTRAN  V 

AUTHOR:  Kim  David  Saunders 

FILE  LOCATIONS:  Absolute,  Reloctable  and  Symbolic 

Elements-  VEGA*L IB .UNORTHOG 

INPUT:  From  terminal  (unit  5) 

Line  1:  NUIN,  NSEG,  NSSEG 
Line  2:  NUOUT 

NUIN  -  Input  unit  No. 

NSEG  ■  No.  of  Segments  desired 
NSSEG*  No.  of  first  segment. 

NUOUT*  Output  unit  No. 


OUTPUT: 


ADDITIONAL 

INFORMATION: 
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APPENDIX  C:  COMPUTER  LISTINGS 


7 

0.  ‘FILE  4  -  TAPE  ID  V002*. 

It  'Profile  iota  logger-Part  2*.ldk  2 
2 s  'LATEST  MOO i  02  DEC  1981 . 164SZIKDS1* . 

3i  'Added  second  lepth  profile  subroutine*! 

4«  *  Added  pressure  Corectlon'i 

5»  ’Changed  profiler  fro*  unit  8  to  unit  5-8  now  res.  for  disk*, 

St 

7|  *0  array  eontrols  output  apt  ions' i 
3»  *  til  eontains  Hazeltine  list  option*! 

9i  '  (2J  eontains  7245A  time  series  plot  option*! 

IOi  *  (Cl  contains  7245A  depth  plot  option*! 

Hi  *  (41  lists  suamaries  of  instrument  motion  statistics*! 

12i  *  (51  produces  a  second  depth  profile  on  unit  703*. 

13t  *  (81  produces  a  plot  of  one  variable  vs  the  other*! 

1A»  'Options  3  and  4  disable  each  other  to  prevent  device  conflict*! 
15i 

16i  'Function  keys*! 

17i  *  f0~  solicit  display  options*! 

18i  *  F0  start/stop  data  logging*! 

19»  *  FI  Generate  end  of  data  file*! 

20i 

21 1  ‘Flags*! 

22i  *  0  -  if  set,  log  to  mag  tapei  if  clear,  don’t*! 

23i  *  1  -  if  set,  close  the  data  file,  if  clear,  continue  logging*! 

24t  *  5  -  if  the  matrices  in  the  DISPLAY  routines  have  been  *i 

25 i  *  -  dimensioned,  fig  5  is  set.  *i 

26 i  *  6  -  if  set,  the  first  pass  through  the  TIMESERIESPLOT  routine*! 
27i  *  7  -  first  pass  flag  for  INSTVEL  *, 

28 i  *  8  -  if  set,  the  2  nd  pass  through  the  TIMESERIESPLOT  routine*! 
29 1 

30 i  di*  A(88, 31,0(9], ZI381 ,T(271 ,Rr0.101 ,S(i:i 
314  dim  Vf 166,201 ,V(10,0,10),Qt(31,U(2],TttB01 
32 i  dim  9(3  31 ,F(4, 31-0(61, EI61 
o3i  'LOAD  VELOCITY  ORTHOGONAL  1 ZAT 1 OH  MATRIX* . 

34.  -.055-»FIl,ll,-.803-»FC2,ll,1.58-*Ft3.11,-.734-»F(4.11 

35.  197-*F(  1  ,21 , 634-»F(2,2 1  ,-.011-*F(3,21, .  759-*F(  4  ,21 

36.  .612-*Fl  1 ,31».349-»F(2,3] ,  -  .  071-^(3,31 , .  309-»F(4 ,31 

37.  trk  Oyldf  5,VI*l,ldf  6,V$,trk  0 

38.  0-»A;ent  'Enter  the  Cast  No.*, All,  11 

39i;ent  'Enter  the  Latitude  in  decimal  degrees* ,AI6, 11 
40 1  ent  'Enter  Longitude-decimal  degrees* , At 7 , 1 1 

41,  cos(S21'*A(  11 , 11 1  -sln(52)-*Al  12,11 

42. 

43i  -3-*0l  Hi  0-*0131 1  0-*0t  21 

44.  on  err  • E RRORHANDLER * 

45.  dev  *nb*,520 

46.  buf  *b2* ,38,1 

47.  buf  ’bl*. 38,1, time  15000, f*d  0 

48.  ’START*. if  f lgl , ldp  7 

49.  buf  *bl* 

50.  cmd  5,*^ _ * 

51.  tfr  *nb* , *bl * ,38 

52.  Jmp  rds(*bl*)"38 

53.  wait  200 

54.  buf  *bt* 

55.  cmd  5  ,*..__* 

55.  tfr  *nb%*b2*,38 
57i  imp  rdst*b£*>a38 

58.  If  flgOygto  *nest* 

59.  cmd  5l*  ,v  * 

60.  cmd  5, *7?*" 

61i  *neat*ifor  1*1  to  38,rdb(*bl*l-*Z(Il»nest  I 
62.  if  Z(ll-79.cli  8,pr t  "ERP'^to  'START* 

63 1 

64,  ell  ’BREAKOUT’ 
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for  !-:  to  66|AtI  ,  1  J-*AI  I  ,21  inaxt  1 
for  I-l  to  38}rdb(*b2*3-*Zl  1 1  .next  ! 

1*  Zlll-79»cli  8ipr t  *ERR*.gto  ■START* 

ell  'BREAKOUT* 

ell  'OPTIONS' 

ell  'DERIVE' 

ell  'DISPLAY' 

if  ab9'.AI15,l]-AI15,31)<ligto  ‘START* 
for  I-l  to  o6}A(l,l]-»AtI,3]'n-xt  1 
g to  'START* 

•BREAKOUT*  « 

■Decode  byte  string  Z  into  prooer  slots  in  A*t 

•PRESS  *  .  :256Zr3Wt213/2  37-At  15,11  .if  bitCO.Zt  1013 |-AC15, ll-AC15.il 

■PRESS  CORR*  ■  .  99897897 ISA [  15 ,11*.  30 325283-* At  15,11 

•TEMP-S*.  <256ZC51*Zt413  /200  0-*At  16, 1 1 

•COND*  t  (256ZI  71 *2161 3 /1000-AC  18,11 

•F-TEMP*  » (256ZI91  -ZIE13/2000-AI17.11 

•SPECIAL  FIX  FOR  GYRE  CPU1SE  ONLY* .Atl7,ll-AI16,l) 

for  1-11  to  27  by  2 1  shf  CZI I  *1 1 ,  -63  ♦  shf  CZI  1 1 ,23-TCIl 

if  bltCO.ZI I 13-1)-Ttll-Tt II 

next  I 

•VELOC  *  «f or  1*1  tc  3}  Tt  U-2C  I -13  1 . 00639-AI22- 1 ,11, next  I 
- AC  20 ,1 1— At  20 ,11 

•I4AGNET*  i for  I-l  to  3»TI  17*2 C 1  -13  1/1000-AC21M  ,  1 1 , next  I 
-At  22, 1 1 -At 22 , 1 1 

■ACCEL  *  .for  I-l  to  3.T123-2CI-13  1/1000-A124*!  ,11  .next  I 
-At 26, II -A 126,11 
•TIME* . 

shf  (Zt  38 1 ,43 -Till  i  Zt  38 1 -shf  (Till  .-43-TI21 
shf  CZ1371 .43-TC31  ,Zl371-shf <Tt31 ,-43-Tr41 
shf  (ZI361 , 43-TI51 ,21361  -shf  CTI5J,-4>-*TI6J 
shf  CZC351 .43-T171  |Zl351-shf <Tt71  ,-4>-TI81 
shf  (Z! 341 ,43— T l 91  }ZC 34 1 -shf  (TC91  ,-4)-Ttl01 
shf  (Zl331,43-TIlll,Zl331-shf  CTt 1 1 1 , -43 -T112 1 
lO0Tfll*lOTt21-TI3]-Tfl31 
10TC41  *TC51-»T  1141 
10T 161  ♦TC71-*TI  151 

10Tt81*TI91*TC10]/13*Tllll/100*TI 121/1000-*Tt 161 

for  I-l  to  4 |TI  I *12 l-»At  1  ♦  1 ,11} next  I 

ret 

•DERIVE*. 


1*2. 

ell 

•SIGMA-T'CAtl6,ll , AC28.11, AC30, 113 

113. 

ell 

'SND  SPEED' 

114. 

ell 

'ANGLE' 

115. 

ell 

'TRANSFORMS' 

116. 

ell 

'TRANSFORM  Q' 

117, 

dl 

'GMAG  CURRENT' 

118, 

ell 

'GMAG  ACCEL' 

1I«  , 

ell 

'GMAG  MAG' 

•27 

ell 

'INST  VEL' 

'  2 '  , 

dl 

'BEST  VEL' 

22. 

if  absCAl 15,ll-All5,31><liret 

23. 

ell 

'dp/dx' 

124, 

ell 

'VEL  SHEARS' 

125, 

ell 

'N  AN0  RI' 

126, 

ret 

•SALINITY*. 

•ARGUMENTS.  1-At 15 , 1 1 ,£-Al 16, 1 1 ,3-At 18, 11  * . 

1 . 45038/9. 9e7-»p4 

6 .  '•658362 1 732e5-p5 .2.0 0529363371e2-p6 
1 . 1 109895 161 2e-2-»p7} -7. 26691983149e-7-p8 
1 . 3586827285e-ll-p9 

p3«(l-5 .25e-6»Cp2-153  ♦plp4)/42.906-»pll 
1 . 60836e-5pl-5. 4845e-10plt2»6. 166e-15plt3-*pl2 


.  .  .  *■  . 


136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 
*55 

56 

157 

158 

159 

160 
161 
162 

163 

164 


pl2/Cl*. 03l786p2*3. 169e-4p2t2)-*pl3 
pll/Cl+pl3)-»pl4 

lC0p2-*pl0»  (p5*p6pl0*p7pl0t2*p8pl0t3*p9ol0t4)(l»tnt(-:6)>-»pl5tpl4/pl5-*pl6 

-  .  08996*28. 8S67pl6*12. 18882pl6t2-  '  '  -61869pl6t3-»pl7 

pf7 *5 . 98624pl6t4  - 1 . 3231 :  pl6t5-*p  1  * 

pl7*pi6(pl6-l) C . 0442p2- . 01 046p2 12- . 004pl6p2)-»pl3 

ret  p!8 


•SIGMA-T*. 

'ARGUMENTS 1 1- AC  16, 1 1,2-4(23.11*1 
•VARIABLES  DESTROYED  >  Ct»I,I,J*. 

If  not  flglOidl*  CC 0i4 , 0 .41 i sf g  10 

for  1-0  to  3 

for  J-0  to  3 

0-*CCI,Jl 

next  J 

next  I 

8.  00969062e-2-*C  C  0 , 0 1 ;  7 . 97018644*- 1-C 1 0  : 1 1 
1.31710842e-4-»CI  0,21 t-6. 11831499e-8-*CI  0,31 

5.8819403e-2-*CC  1 ,01  j  -3. 25310441e-3-*Ct  1 ,11  |2.8797153e-6-*CCl  ,21 
-8.H465413e-3-*CI2 ,01  .3.89187483e-5-»C(2 ,11 
4.  7£600414e-5-*CC3,01 
0-*p3 

for  I -  0  to  3 
for  J-0  to  3 

if  I*J<4ip3«C(I ,  JI*oltI«p2tJ-*p3 
next  J 
next  I 


165 

166 

167 

168 

169 

170 

171 

72 

73 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 
165 
186 

187 

188 

189 

190 

191 

192 

193 

194 

195 
136 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 


ret 

•SND  SPEED* > 

•SOUND  SPEED  FORMULA* i 

AI 16 , 1  J-*X  j  A 1 28, 1 1-»S  i  Al  15 , 1 1-*D 

100C 1449*4. 6X-.055XX*.0003Xt3*(l. 39-. 012X1 CS-35) ♦. 017D)-*C 

C -*4(29,11 

ret 

•TRANSFORMS*.  At  49,  U-*pl 
deg  jcoaCpl  J-*p2j -alnCpl)-*p3 
0-»X-*Y 

for  1-1  to  3 1 At  1*21 , 1  lt2*X-*X| Al  f  *24, llt2*Y-*Y 
next  1 
FX-*X  ,-TY-*V 

for  1-1  to  3|A(  I  *21 , 11  /X-*EC  1 1 1  AC  I  *24,1 1/Y-*E  1 1*31 
next  1 

for  1-1  to  3|-E(  I*31-*B(3, 1 1  inext  I 

for  I-l  to  3 1  (EC  I  l-p3B(3, 1 1 1  /p2-*B(2 , 1 1  inext  I 

xgn(A(22, 11 JFC1-BC2, 1172-8(3,11723  ■♦8(1,11 

A(26 , 1 1 A(24 ,1 1 -A( 27 , 1 1 AC23, 1 l*pl 

-B(l, 1 1 (A(25, 114(24,11-4(27, llAC22,ll)/pl-»B(l, 21 

Bt 1, 11 CAI25 ,114(23, 11-AC26, 11AC22, ll)/pl-»B( 1,31 

ret 

•TRANSFORM  8* » 

for  J-l  to  3jF[  1 ,  J1-*X 

for  1-1  to  3tACI*18,llFCI*l ,  J1*X-*X ;  next  I 

X-»A(  J*30 , 1 1  j  next  J 

ret 

•GMAG  CURRENT*. 

•CURRENTS  IN  GEO MAG.  COORD.  FROM  INSTRUMENT  COORD. •. 

for  1-1  to  3 

0-*Y 

for  J-l  to  3|BC1 , JJAI J*30,ll*Y-*Yinext  J 
Y-*A(  I  *39,11 
next  I 
ret 


GMAG  ACCEL* . 

ACCELERATIONS,  LESS  0,  IN  GEOMAGNETIC  COORDINATES* • 


»07.  for  1*1  to  3 1 Q-»X 

208.  for  J*1  to  3<Bt  I ,  JIAt  J*24 ,1  l*X*X|ne*t  J 

209.  X-*4C  1*33,11 1 next  1 

210.  A 1 36 , 1 1 *9 . 9-*A (36,11 

211.  ret 

212. 

213:  ‘GMAG  MAG*. 

214.  ‘MAGNETIC  FIELD  IN  GEOMAGNETIC  COORDINATES* . 

215i  for  1*1  to  3j0-*Y 

£16 .  for  J-l  to  3|BC  1  ,J3AtJ*21 ,  ll»Y-'Yjne*t  J 

217.  Y-*AC  I *36,1 1  .next  I 

218.  ret 

219. 

220.  ‘ANGLE*. 

221.  ‘COMPUTATION  OF  ANGLE  BET.  G'  AND  H»*. 

222.  TCAt22,llt2*A[23,llt2*At24,llt25-»H-*AC5Q,ll 

223.  A(22 , 1  J  /H-»0I 11.4(23,1  J  /H-*DI2  J  » At  24 , 1  J/N-*0I31 

224.  JCAI25,llt2*AI26,llt2*4t27,llt2>-HS-*Al51,il 

225.  AI25,ll/G-»Dl41,AI26.11/G-*DC51|Al27,ll/G-»Dt61 

226.  0-»P|for  J*  1  to  3iD( J1DC J*31*P-*Bjnext  J 

227.  oanCB) -»B  jB-»Al49, 1 1-*A19 , 1 1 

228.  ret 

229. 

230.  ‘INST  VEL‘« 

231.  ‘INSTRUMENT  VELOCITY*. 

232.  AI2 , 1  l-»pl 

233.  (AI 2 , 1 1 -pi  186400*360 0A 1 3: 1 1  *60At  4 ,11*4(5,1  ]-*p3 

234.  CAI  2 , 21  -pi  186400*360041 3, 21  *60 Af  4.21  *AI5 , 2  l-»p2 

235.  if  f Ig7,p3~p2-Al55,lj 

236.  if  At 55 , 1 1  -  0 i5-»AI55 , 1 1 

237.  if  not  f  lg7  |5-*A(55 , 11  j  sfg  7 

238.  l00CAri5.11-Atl5,21?  /  A  155 , 1 1  -*At  45 , 1 1 

239.  0--AI 43, 1  l-*At 44 ,11 

240.  ret 

241. 

242.  ‘BEST  VEL ■ . 

243.  ‘CORRECTS  OBSERVED  VELOCITIES  FOR  INSTRUMENT  MOTION*. 

244.  for  1-1  to  3<  At  I  *39, 1  J-At  I  *42 , 11  ■♦At  1*45,11  inext  1 

245.  ret 
246> 

247.  ‘dp/dz*. 

243.  ‘VERTICAL  DENSITY  GRADIENT  Cdp/ds>*. 

249.  At  15,  II -At  15.31-.AI56.il 

250.  (  At  30 , 1  ]  -  A130 , 31 J  /  AI56, 1  l-»Af  57, 1  J 

251.  ret 

252. 

253.  ‘VEL  SHEARS*. 

254.  ‘VERTICAL  VELOCITY  SHEARS  *. 

255.  for  1-46  to  48 1  C  A 1 1 , 1 1  -  A 1 1  , 31 5 /At  56 , 1 1 -»At  I  ♦  12 , 1 1 

256.  ret 
157. 

258.  *N  AND  RI*. 

259.  *N  AND  RICHAROSON  NUMBER*. 

260.  980At57.il/Cl*  At30 , 1 1 /l  0001-.pl  i  if  pl<0,-pl-*pl 

261.  ■fol-*AI62 ,11 

262.  A 1 58, 1  lt2*AI59, 1  It2-»p2 

263.  if  p2#0ipl/p2-*AI61 , 11 

264.  if  p2-0 j999-»A161 , 1 1 

265.  ret 

266. 

267.  ‘DISPLAY*. 

268.  ell  'HAZELTIME* 

269.  ell  ‘TIMESERIES* 

270.  c 1 1  'DEPTH* 

271.  ell  'DEPTHS* 

272:  ell  *V1V2* 

273.  ell  ‘STATS* 

274.  ret 

275. 

276.  ‘HAZELTIHE*  .if  Ot  1 K0  i abeCOI  1 1>-*0fll  ,cl  1  'HAZIHITI AL IZE* 

277.  if  Ot 1 1-0 i ret 
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278. 

279. 

280. 
281. 
282: 

283. 

284. 
285  . 
286. 
287. 

285. 

289. 

290. 

291. 

292. 
293  > 

294. 

295. 

296. 

297. 
299. 

299, 

300, 

301, 

302, 

303, 

304, 

305, 

306, 

307, 

308, 

309, 

310, 

311, 

312, 

313, 

314, 

315, 

316, 

317, 
313, 

319, 

320, 

321, 

322, 

323, 

324, 

325, 

326, 

327, 

328, 

329, 

330, 

331, 

332, 

333, 

334, 

335, 

336, 
337: 

338, 

339, 

340, 
'.41, 

342, 

343, 

344, 

345, 

346, 

347, 

348, 


ell  'HAZWRITE' 
ret 

•HAZINITIALIZE*. 

if  not  flg5,dim  Ltt3,801,afg  5 

•  ‘’imi 

for  1-1  to  VCOtll, 01, VttVCOtll, II, 1,61+LOtl, 1-6(1-13,611, next  1 

wtb  3,27 ,26, wtb  3,27,17,0,3 

fmt  l,c80,zjwrt  3,LtCll,ret 

■HAZWRITE* , 

fmt  1 , f 6. 2 , z ,  fmt  2,/ 

wtb  3, 27, 17,0, 5, wtb  3,27,26 

for  1*1  to  Vt0tl],01,wrt  3.1, AC VI Oil), II, 11, next  1 

wrt  3. 2 .wtb  3,27,17,40,2 

fmt  3,f3.0,2x,f2.0,2x,f2.0,2x,f6.3 

wrt  3. 3, AI2.11, ACS, 11, At 4, 11 ,AI5,11 

wtb  3,27,17,2,2 

fmt  4 , "CAST  •  *  ,f 4.0 jwrt  3. 4, All, 11 

wtb  3,27,17,2,l»fmt  5, * LATITUDE  *,f6.3,*  LONGITUDE  *,f8.3 

wrt  3. 5, AC6.11, AC7.11 

ret 

•DEPTH'. if  OC31<0,-Ot31-QC31, 0*0141, ell  'DEPTHINITJ ALIZE* 
if  0C31  **0  i  ret 
ell  'DEPTHPL0T* 
ret 

'DEPTHINITCALIZE'.wtb  706,27,85 

wrt  705, 'IP, 1000, 1000, 6000, 6000'ipac  705,pclr 

f  xd  0 

ca i z  3. 5 , sc  1  0,10, 0,10, pit  0,16,l,lbl  T* 

caiz  2,acl  0 , 1 ,300 , 0 , xax  Ojxax  303, yax  l,20,yax  0,20,0,300,5 

caiz  3, 2, 1,90, pit  - . 1 , 170 , 1 , lbl  'Pressure* , cai z  3, 2, 1,0 

for  1-1  to  VCOC31.01 

/a  1  (V$  C  VCQC31 ,11 , 7, 1313*pl 

valCVt  CVC0C31  ,  I  1 , 14 , 20 1 1  *p2 

AC  15,1 1 *R 1 0  1  ,ACV1QC31 , 1 1 ,11  *R 111 

5-»p6,if  abs <p2-pl) <10 , l*p6 

l-»p7«if  absCp2-pl)>99,50-*p7,  l-*p6 

caiz  2 , ac I  pi ,p2 ,0 , 10 , xax  10* 1 ,p7 , pi , p2 ,p6 

caiz  3, ac 1  0,80, 0,10, pit  30,10.1-1,1,161  VtlVtOCSl , 11 ,1 ,61 

caiz  1.5 

next  1 

ret 

■DEPTHPLOT'.pac  705 

for  1-1  to  VC0C31 ,01 

val<Vt(VC0(31,I],7,1313*pl 

val CVt(VIOI31 ,11 , 14,201>*p2 

acl  pi , p2, 300,0, 11m  pi, p2 , 300,0 

ACVC0C31,!]  , 1 l*p3, At  15, 1 l*p4 

pit  p3,p4 , 1 ,pl t  p3 ,p4 , 2 

p3*R  til,  p4*R  C  0 1 

next  I 

1  im 

ret 

•STATS*  » if  0141-0, ret 

if  Ot  41  <0 ,  -0t41*0t41 , 0-»0t31  ,cl  1  *  STATS  INITIALIZE* 
ell  *STATSXEQ* 
r  et 

•STATS INITIAL IZE* , wtb  706,27.85 
wrt  706, T« 

for  I-l  to  7 , 0-»Sl  1 1 ,  next  I 

fmt  1 , 20x, *  INSTRUMENT  ATTITUDE  STATISTICS* , / 

fmt  2,5x, "PITCH* ,14x, •VAN* , 13x, *NAO  DIP* , 13x, •PRESSURE* 

fmt  3,*  Mean  Std  Dev  Naan  Std  Dev  Naan  Std  Dav  Naan 

wrt  706. 1 .wrt  706. 2, wrt  706.3 

rat 

•STATSXE8* iSCll-l*St  11 

•Pitch* ,  dag  lacaCBt  3,31  l-*X,e  11  *SUPDATE1*C2,X> 

•Vaw*  .acaCBtl  ,lll-*Xtcll  *SUPDATE1*C4,X) 

•MagDip* , At 49, 1 l*X,c l 1  *SUPDATE1*C6,X1 
•Praaa  •  . AC  19, 1 1*X ,e  1 1  'SUPDATE1  * C8.X1 


if  sm-2«f9te  *SLIST* 

*SLIST*.for  1-2  to  8  by  8190b  •SUPDATE2* 

BO*t  I  ’ 

fat  l,f7.2,2x,Z]f»t  2,/ 

for  1-2  to  8f*rt  706. 1 ,Sll I i»oot  Ijwrt  706 

for  1-1  to  7(0-*SIIliBOot  I 

rat 

■SUPDATE1*  iSIoll«p2-»Stpt  1  |S(pl*n«p2t2-»SCpl«l  1  trot 
•SUPDATE2*  tST  1 1/SI11-*SI  I J  |<TC£{  1  *1) /St  11 -St  1  lf2)*SI 1+11 trat 

•TIMESER1ES*  1  if  01 21<0» -0121-»0C21 1 0-*OI51  tell  'TIMESERIESINITIALIZE* 

if  0121-0 • rot 

ell  'TIMESERIESPLOT' 

ret 

■TIMESERIES INITIALIZE* iwtb  704,27.85 
fxd  Oicfg  8 
pi:  703ipclr 

art  703, MP,  1000, 1000, 7500, 10000* 
csiz  2 

sc 1  0,2000,0, 1| too  0,100,0, 2000, 5iyax  2000ixax  Ijyox  0 
VCQI21 ,01**pl 

for  I-l  to  pl-ltxax  I/pl , 100 inext  I 
csiz  2,2,1,90 

for  I-l  to  pltplt  -200,CI-l>/pl*.02,l»lbl  V0IV1OC21 ,1 1 , 1 ,61 inext  I 
9000/pl-»p2 

for  I-l  to  pljfat  l,*IP*,*,*,fz4.0,*,*,fz4.0,*,*,fz5.0,*,*,fz5.0 

wrt  703. 1,1000, 10 30*(I-l)p2, 7500, 1000+Ip2 

csiz  pi, 2, 1/6. 5,0 

val  (VO  1 VI 0(21 ,11,7, 131)-*p3 

v«l(V*rvt0t2J,l )  ,14,20))-»p4 

9 c  1  0 ,2000 ,p3,p4|f xd  0j5-*p6|if  sbs(p4-p3X10|l-»p6 

l-»p7(lf  obsCp4-p3}>50t2->pb|50-*p7 

if  Iaod2 iyox  0 ,p7 ,p3,p4 ,pC 

if  not  I«od2 iyox  -80 , p7 ,p3 ,p4 ,p6 

next  I 

ret 

■TIMESERIESPLOT* ipse  703 

if  not  flo6  or  flg8>«to  ‘TSlnext* 

AI3,ll->Ql3]tAI4,ll-QC4]|AC5,ll-»QC5I 

A 1 2, 1 1-*QI  1 1  ;3600AI3, 1  l*60At  4,  1]*AI5,1J-*Q[2J  j  sf  g  8 

wrt  703, *IP, 1000, 1000, 7500, 10000'tscl  0,80,0,10 

fxd  Otcsiz  1.5 

pit  10,10.3,l|Ibl  'Start  Time  ',013],*  .*,8141,*  .*,0151 
pit  10,10.5,1] lbl  'Start  Day  *,0111 
■TSlnext* 1 

if  flgGigto  *TSnext' 

dim  QI71|AI2,1 1-*QI  1 1  $3600 A 13 , 1 1+60AI4 , 1  I  *At5 , 1 1-*®I21 
sfo  6 
■TSnext*  t 

0-*p8»if  AI2, 1 1>0I1 1  j 86-40 0-»p8 

V10I21 ,0I-*pl|90C0/pl-»p2 

p8+360  0AC3,1 1  *60  AI4, 1J+AI5, 11-012  J-»p9 

for  I-l  to  pljfmt  l,*IP*,*,*,fz4.0,*,*,fz4.0,*,*,fz5.0,*,*,ft5.0 

wrt  703. 1,1000 ,1000* (I -l)p2, 7500, 1000* Ip2 

val  (Vf  FVI0I21 , 1 1 ,7, 1317-*p3 

val (V0 1 VI 0121 , II  ,14 ,201>-»p4 

9c 1  0 ,2000 ,p3,p4| 1 in  0,2000,p3,p4 

Atvtoi2i,n,ii-*pie 

pit  p9,pl0,i 

pit  p9,pl0,2 

next  I 

11m 

rat 

-ERR0RHANDLER* 1 if  arn-4,prt  'TIMEOUT  ERROR* 

prt  *ERN  -* ,arn 

prt  *Lina*,arl 

prt  *R0M*,roa 

tine  0 

cli  7 

ell  5 


420 

421 

422 

423 

424 
125 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 
159 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 


gto  •START* 

•OPTIONS* i l f  A-0,ret 

cnt  'Enter  Display  Title*, TO 

•YESYES'.ent  'Enter  Display  Device  No.*,A|if  A-0,ret 
If  A)9  or  A<0 | jap  -1 

ent  ‘Enter  display  option  no.*,0CAl|lf  Ot  A]>0?-OIA)-*OIA1 
ent  *Do  you  want  wore?' ,Q$ , if  cepCQOt 1 , 1 l)-*Y* ,gto  "YESYES* 
;f  capCQOll , 11)#*Y*  and  capCQOI 1 , 1 1 )e*N* , jap  -1 
0-»A,ret 

•SIGTEST'  .fat  l,3fl0.5 

ent  'enter  T',T,ent  'enter  S',S 

ell  *SlGHA-T*CT,S,a) 

wrt  3. 1 ,T,S,Q 

gto  'SIGTEST* 

•TSTST'iell  'T1MESERIESPLCT' 
stp 

•DEPTN2*  t If  OC 51  <0; -0151 -*015]  ,0-»0121  jell  'DEPTH1HITI ALI2E2* 
If  0151*0, ret 
ell  '0EPTHPL0T2 ' 
ret 

•DEPTHIHITIALIZE2'  >wtb  704,27,85 

wrt  703, 'IP, 1000, 1000, 6000, 6000' ,psc  703,pclr 

f  xd  0 

cslz  3.5,scl  0 , 10 , 0 , 10 , pi t  0,16,lilbl  TO 

csiz  2,scl  0 ,1 ,300 , 0 i xax  Ojxax  300iyax  1,20, yax  0,20,0,300,5 

cslz  3, 2, 1,90, pit  - . 1 ,170 , 1 , lbl  'Pressure' , calx  3, 2, 1,0 

for  1*1  to  VI0151 ,01 

val  (V0 1 VI 0151 ,11 ,7,131)-*pl 

val  CVttVtOISl  ,Il,14,2011-*p2 

Atl5,ll-»Rf01,A[Vt0[5],Il,ll-»R[l] 

5-»p6,  if  obs<p2-pl  J  dO  ,  l-*p6 

1"*P7  ,lf  ebs(p2-pl  1  >99,50-*p7 1 l-*p6 

csiz  2, scl  pi ,p2 ,0 , 10 i xix  10+1  ,p7,pl,p2,p6 

csiz  3, sol  0,80, 0,10, pit  30 , 10. 1* i , 1 , lbl  VtlVtOtSl , 1 1 ,1  61 

csiz  l.S 

next  I 

ret 

'DEPTHPLOT2* apse  703 

for  1-1  to  V10I51 ,01 

val(V01V(OI51,U,7,131)-»pl 

val CV01VI0151 , 1 1,14,20 l)-»p2 

scl  p! , p2, 300,0, 1 lm  pi, p2, 300,0 

fitVI0l51.il  ,l?-»p3, At  15,ll-*p4 

pit  p3,p4.t,plt  p3 ,p4 , 2 

p3-»Rl  1 1  ,p4-»Rt01 

next  1 

1  la 

ret 

•VlV2*«if  Ot61<0,-Ot61-*Ql61;cll  '7IV2INIT* 
if  CI61-0,ret 
ell  'V1V2PL0T* 
ret 

'V1V21NIT*  »ent  'Enter  the  plot  unit'.Z 
if  Z*703,0-*0t51-*0121 
if  Z*705.0-*0l41-»0t31 
wtb  Z*1 ,27 ,85 
fxd  Oipsc  Z 

wrt  Z,* IP, 1000, 1000, 6000, 6000* 
ent  'Variable  number  l'.UHl 
ent  'Variable  number  2*,U(21 
csiz  2 


484.  val  CVttUt  11 ,7 , 1317-.pl  ,ve  1  (V$(0(  1 1 ,14 , 2311-*p2 

485.  val  CV*  101 21 ,7 , 131  >-»p3,  val  CVttUI  21 ,14 ,20Il-*p4 

486.  scl  pi  ,p2,p3  ,p4,5-*p6,  if  abs(p2-pl><10,  l-»p6 

487.  l-»p7,if  ebs<p2-pl)>99»50-*p7,l-*p6 
438.  xax  p3,p7, pi ,p2 , p6,xax  p4 

489.  5-»p6, If  ebs tp4-p3> <10 , l-»pC 

490.  l-*p7,if  nbs<p4-p3>>99,50-*p7,  l-*p6 


491 «  yax  pl,p7,p3,p4,p6,yax  p2 

492 i  caiz  3, ad  0,80, 0,10 iplt  30 , 11 . 1 , 1 , lbl  V*IUI1J,1,6) 

493 ■  c 9 1 z  3.5|9cl  0 ,10, 0 , 10, pit  a, 15,1, lbl  T* 

494:  9 c i  0,10, 0,80, cals  3, 2, 1,90, pit  -1,30,1, lbl  V»[Ut2],l,61 
495 i  C9lz  ,lim 
496:  ret 

497.  'V1V2PL0T* :p9C  2 

498:  val  CVtCUIll  ,7 , 1311-»pl ,  val  (VtIUI  11 , 14 ,20  J7-*p2 

499.  val  CV$tUI21,7,13n-*p3,  val  CVS  (U121 ,14,20  l>-*p4 

500.  9 c 1  pi ,p2,p3,p4 • 1 iw  pl,p2,p3,p4 

501.  AtUlll.ll-pSjAtUm  ,u-»Pe 

502.  pit  p5,p6,l,plt  p5,p6,2,lim 

503.  ret 
■9743 


IN  DISPLAY  OPTIONS 

OPTION 

/  1 

OPTION  t 

1 

OAY 

1  PRESS 

o 

HP 

2  Ul-GC 

3 

MIN 

3  U2-GC 

4 

SEC 

4  U3-GC 

5 

CAST 

5  ABSG-0 

6 

LAT 

6  N 

7 

LONG 

7  R I 

8  T-SLOU 

9  SAL 

OPTION 

*  2 

1 

PRESS 

2 

Ul-GC 

OPTION  • 

3 

U2-GC 

1  T-SLOW 

4 

U3-GC 

2  SAL 

5 

UI3-GC 

3  SIGMAT 

6 

ABSG-0 

7 

ABSH-0 

8 

MAGDIP 

OPTION  # 

1  Gl-IC 

2  G2-IC 

OPTION 

•  3 

3  G3-IC 

1 

r-SLOH 

4  Hl-IC 

2 

PRESS 

5  H2-IC 

3 

Al-GC 

6  H3-IC 

4 

A2-GC 

5 

A3-GC 

6 

Hl-GC 

OPTION  # 

7 

H2-GC 

1  Ul-GC 

8 

H3-3C 

2  l'2-C-C 

OPTION 

*  4 

1 

PRESS 

2 

Ul-lC 

3 

U2-IC 

4 

U3-1C 

5 

Ql-AP 

6 

G2-AP 

7 

Q3-AP 

9 

THETA 

9 

PHI 
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iAK»*u ‘Transcribe - _ _ 

l  c***************** ****************************************  ********* 

-2 - C*************  *********  *********** ****************.******** »»»«»+»»» 

3  C***  PR CGR AH  t  TRANSCRIBE  *** 

r  »- transU-ifs  pr  of  her  data— to. .engimeep 

5  C***  LEAVING  VELOCITY  DATA  IN  NON-OR TMOGONAL  *** 

6  - C*** - INS-TRUN£HT-  COORDINATES. - *** 

7  C***  #** 

* - C***  AUTHOR — X— J.J •— VEG  A  ,-COMPUTER SCIENCEGCORPv-CPRIHAR  Y» - *** 

9  C***  K.D.  SAULOERS  .  NORDA  (  SECONDARY  -  REQUIRED  IN  *** 


DOCUMENTATION.) 

DATE  IOF  LATEST  REVISION)  t  )5  MARCH  X9S2 


INPUT  s  FILE 


TYPE  OF  OATA/COMMEnTS 

I  IF  USED  FRO**  A  TERMINAL,  THE  PROGRAM 
— SUPPLIES- SOLICITATION  PROMPTS  » . — 

data  tape  in  hp  9«2Sa  /  norda  format 

- L--T-HE — T  APC  SHOULD— BE  COPIED-FROM  THE— 

BOO  BPlt  FAW  TAPE  TO  A  HIGHET  DENSITY 


PROGRAM  WILL  RUN  FASTER  Af;0,  second, 
-YDU-WILL-  NOT  GET-ARNORHAL-FRAME  COUNTS- 
WHICH  WILL  TERMINATE  THE  PROGRAM?  THE 

-PRO  RED— METHOD— E-OR— COR  ylNS—XS— USE - 

RCOPY.MN  INPUT* 10  .) 


C***  OUTPUT  *** 

-C*  ** - *** 

C***  20  FEB  FILE  CONTAINING  PROFILER  DATA.  *** 

C*** - l-TH-IS— FILF-  MUST  GE-ASSIGNEO-PRIOR  TO +** 

c***  executing  the  program.)  *** 

-C*-»* _ ***_ 

C* ************************ ***********v** ******************  ********* 

c***** ********  ****************************************************** 

c 

C - 

C  MAIN  CALLING  ROUTINE  FOR  TRANSCRIBE  ROUTINE 

_C _ 

c  this  program  reads  oata  from  mag  tape, 

C  DECODES- IT,  AND-STORES  -IT-IN-  FEE  FILE-S* - 

C 

-C - 

C  ARRAY  LIST 

-C - OBLK - IS— J.HE— INPU-T— DAT  A  ARRAY _ 

C  READ  FROM  TAPE. 

C - RD - IS  THE-  ARRAY  CONTAINING _ 

C  THE  BYTES  OF  THE  OATA 

C - STRTNg  TO  BE  DEOCOEO. - 

C  VAR  IS  THE  ARRAY  CONTAINING 

_C - IHE-Of  CODED -vArLAELES  OE _ _ 

C  A  data  STRING. 

C  - VW - IS  THE  FEB  TILE  OATA.  ARRAY _ 

C  ACOC •  IS  THE  FEB  FILE  ALPHA- 

C  - NUMERIC-HEADER  ARRAY-. - 

-C - FPOcN  -IS— T-Hg— F£6-fiL-E— FL-DA TING - 1 - 

C  POINT  HEADER  ARRAY. 

-C - I-OOCN— IS — THE—  FES-FTL-C-INFF-GEr - 

c  HEAOER  array. 


C  SEE  KAVO  tech  NOTE  'THE  FAST  AND 


23 


-€A*-7— Bt-NARV-IFEB-I  F  ILE1  R* - 

2.R .  HALLOCK  FOR  MORE  ON  FEB 
FIt-tS-. - 


•nrrrimj 


OONMON/NHOR /LU , NNrNBN ,N* 
COMMON/W0ATa/VW( IS,530> 
COM  MON/fcOOC  f/F00C*<40-4— 
COMMON/ WOOC I/IDoCW( 20  I 


•?«]'!’**» ’fj*!  fAflfi.H" 


-C - 1NIT4AL-«E-CONTrOL-HCAO£R - 

C 

DATA  tlpWII »,I=l ,15)/ •PRESS*.,STEKP’,*CONDt ,»FTEM1»» , 

- l_-t¥LQcU-V,VLaci2Jr1-w*-0CK?^*UPW-i,-,-*A64^-,NA£«i»-!ACCLli-** 

S  ’AC CL  12  ’,  ’ACCLI 3  *,»  TIKE  •  ,  'RELSEC  •  / 

- -CA-'PA— AiNNWyNNIP-fNAlF-t-NN-Jt  l«NA/FOO^15,44)r20,  SO/- - 

C 

r  on  n  inr'p  inoiw 


NR  IT El  6,  5030  > 

- read — <s-,ior  iaoo-cwi ia,i-i-»74 - 

Read  is.ioi  i aoocniii  ,i=i3f iri 

- READ  -  <S  »1Q)  C  ADO  CM  (-1 }  25-,-I-U - 

C 

C  READ  ID  DC  ARRAY 

KRITEI6, SDDl I 

- read — is,40i  -ioocauza.,  ioocr4_U-,jdocu  ipA^ioocui  iot,idocu(  iaj — 

c 

-.C—AfEAD-FO  oC-AftRAy - - - 

writei 6, 5002  * 

WRITflfc, 50031 

- READ-AStAC)  ( FOoCh(  1)  ,ls|«v4*| - 

WRITE»6.5004» 

- REaD_45.40»  IrDOCWC-U  ,1-21.-274 - 

WRlTE<6.5005» 

- READ  U,4Gl  IfCQOUU, 1=31, API - 

C 

C  RE  AO  FILE  NAME,  SEG  RENT- NAME  , -HSGu-LEV  EU, - 

WRITE (6, 5006) 

- J:E*D  IS.5G>  HMFM - 

W'R  IT  El  6, 5C07 1 

- lead  is,s;.)-n*£n - 

WRIT  El  t* 5008  1 

READ  IS  ,40)  HSGw - 


-FORMaT-ITAOI 
FORMAT  II 
-FOR  MA-T — IAO-)— 


A00C«138»=*  COORO* 

- ADOCV  (59)  =  ’ - *- 

CALL  SUBROUTINE  TRAN  SC 
- CALL-TRANiC - 


FORMAT! ♦  INTER  3  LINES  OF  ALPHAMERIC  DOCUMENTATION*! 


4 


126  - SOW - FCRMAltJ~£NTER-AG-lNTCGeRSACPttl$£-N0.,S7  fit  NO. ,*r~ 

127  1  ’ABSOLUTE  NO.  Of  1ST  SEG.tlNPUT  TAPE  NO.  .YEAR  of  CAST*» 

128  - 5(10* - F  ORNA-TT-i— ENTER- (FL0-ATIN6--PT.-I — *-*•/- - - 

129  1  •  SAMPLE  INTERVAL  ISEC),  2EP0»  ST*RTI*S  :ATjTUOE  COEC.  DEG.l* 

_130 _ 2 — t— S-TAJtT-lNft— LO^GT  TUDE-t— T-I-^  OF  FIX  ■  AT— M-AR-T — (DEC-.— 0A-YC4-1/ - 

131  3  ‘  ENDING  LATITUDE,  ENDING  LONGITUDE,  eNqINg  T iNf ,  •/ 

-  132 - A — *- MaXImUm-pRESSUpE— ID£CIBArs>— *-*/* —  . - —  - 

133  5003  FORMAT!  •  ENTeP  (F»  S  MAGNETIC  VARIATION  »ND  DIP,*/ 

- ..  1 3# _ X— ?-S  HI  P-  SPEE  0  AND  HE  AO  IN  G -3  /  /  > - - - 

135  5009  FORMAT!  •  ENTER  DRY  BUlB  TEMP.,  NET  B»L  '  TEMP. , SURFACE  TEMP. 

__14G _ 1  /*  BABOMETR  IC — PRESSURE,  -W  JN-p-SPEEQ  -4N-TS4-. .  N^AEO— 0 1  HOC  T  lOt,.  ■>  / — 

137  2  *  SIGNIFICANT  W»«E  HEIGHT  (FT»  »//> 

138  _ SOOS _ F-OCEHaTL’- ENTER— (F-J— :-CASI_St*R-T-DAY,  HOUR^-MlH»$CC*/ - 

139  1  •  ZERO,  CAST  END  TIME:  DAY,  HOUR,  HIN,  SEC* //I 

1H0 _ SOOG _ F-ORMATt-t-ENlER  -THE-sEOUENTlAL— FILE  NUMBER-ID  •  I - 

HI  S0U7  FORMAT!  •  ENTER  NMe«  WHICH  CONSISTS  OF  A  3  DI3IT  CAST  N0»»* 

_ 142 _ 1  *  AlO  A  1  DIBIT  PaQ£JLE_Nfl-_!_^» - 

193  5008  FORMAT! •  ENTER  THE  MESSAGE  LEVEL  C-9*//» 

149  __ - END - 


/ 


« 


RPRT.S  V.TRANSC 


1A»LIBU  )  ►TRANSC- 
1 

- 2- - C - 


SUBROUTINE  TRANSC 


3  C  THIS  ROUTINE  READS 


s  c  oecoce  each  data 

* - c — CYCLES-AR  e- decoded, 

7 

8 
9 

-w- 


S INGLE  DATA'  BLOCKS  FROM 
5 U»R OUT  IN6-5 — TO  FIND  ANP- 
STRlNG.  AFTER  500  DATA 


-F-E6— SEGMENT— I -S — wR-I-TTtN-»- 


- — C  — COM  HONS-  P0R-2WR-W- 

c 


-Com  ho  n/  whom  imm******  .nhbw,  mhfn  ,-wfw  .newsman  ,»  pn  n  «  i 


n 

12 

13 

-IR¬ 

IS 

16- 

17 

IB 

19 

20- 

21 

-22- 


COMMON/WOATa/VW! 15,57QI 
-COM  MO  N/ vOOC  F/ FOO  Cfc  1 RO-1— 
COMMON/ iiOOC I/ID0C*I2C  ) 

— COMMON/ wDOC-A/aOOcH!SOT — 


-COM HON/WA&WHS b P , H  5G  W , NNNR  ,NNHW  ,NN IP  rNNf ,N  NT ,N  N  A , T  RST rtWST- 


COmhON/EETA /BHAT !3,3) ,0RTH0GI H , 3 » ,A CCAL ! 3 > 

-INI  T1AL17E  VELOCITY  ORThOgONALIZATTON-H^TRIX- - 

CAT  A  OR  TH0G 1 1 , 1 1 ,0R THOG |2 , 1) , ORTHOG 13 , 1 1 ,0R THOG ! 9 ,1 )/ 

- B-~,055  f-»003-»l  «58  ,-.799/- - 

DATA  ORTHOG (1,21 .ORTHOG (2 ,2 > , ORThOG (3,2 ) .ORTHOS (9 ,2 » / 
4-TT,--.A8,  »-»0Hr«~T3TA  ■ 
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-29-  - 
25 
26-- 
27 

-26— 


data  ORThOGI1,3).ORThOGC2,3»,ORTH06(3,3I,ORTHOG(9,3I/ 
-,612  ,  ,8-99,  —  ,&7 1,. 8 09 /  1  —  ■  —  * 


— -INTEGER— 8L-KSI 7- 


OATa  BLkS!2/99S/^ 


29 

30 

31 

32 

33 
39- 


OlMENSI ON  VAR (17  I*Rd! 3TI 

- 1NTE6ER-0BLK I  500  » - 

DOUBLE  PRECISION  OP  TIME .RELOAY 
- RE  AL-W (31, LENGTH - 


kflg=i 

-lOOCWIM-O- 


35 

36 

37 

38 


IOOCn  |4 )=0 
-NHBYtE=1- 
N2-0 

-IHRSTr  1- 


25 


39 


KflAD-Q 


90- 

91 

92 

93 
99 
95 

-  A6- 

97 

98 

99 

50 

51 

52 

53 
S9 

55 

56 
ST- 

58 

59 

60 
61— 
62 


_l5EC-=-0- 


C  - 

c 

c 

10 

-C— 


NBDt0I  =  0 


*9*9* 


CONTINUE 


IF  (KFLG.EQ.Ol  GO  TO  30 
C  - 

C  99999 

c  rcao  next  input- slock^ - 

20  continue 


_410  - 1  GOO — L-  -  1-  ,600 


li'CO 


CPU  (II 
CONTINUE  — 

kflg  -  o 

CALL  NT  RAN  -  (IQ  •  2  9»  OBLK  tL  l 


-IF— (l  .CO.— W~ CAL-L-NTRAN— t-iO-r92 I- 


IF  iL.E0.-2)  60  TO  22 
-IF — IL-.E  Q.-3.0R.L-.E0.-A h  SO.-TO  1020- 


KFL6=0 


60  TO  30 


999  99 

IF-X  CE_fcE  AO.  .-WRITE  PARTI  A  L-SEGHEN I . A  NQ-JgJmi  AT£~RO  U  T  JN  E, 

NFLG=1 

- 100CM  U  1  =  1-  - - ; _ 


6A 

65- 

C 

c 

66 

--67— 

22 

68 

70 

71 - 

- C— 

72 

C 

-73 _ 

_ C 

WRITE  (6,23  1 

TOR  MAT  U-EOE-OEAP-  BEEP RE-SPECIFIED  £N0  OF—CASt-TlM 


*  *  *♦  PARTIAL  SEGMENT  WRITTEN  99*1 


79 
75- 

76 

77 

78 
-79 

80 
81 
«2 
83 
89 

-85- 

86 

87 

88 
89 
93 
91- 
92 
93" 
99 

95 

96 
97- 

98 

99 
100 


99999 

Mil  boiiiimf  TO.  read  Ann  ncrnnr 


30  CONTINUE 

ALL-OA  TSTR — 4 XFL- G,0 BLN, rLkS  I2-r8NB-TT-£rR0 ,  VAR  I- 
_IF — (KFLG.EO  .L — »  AND.  NFL6 — .NC^-14 — GO— TO— 20 - 


C 

-C- 


-99*99. 


c 

c 


C  OECOOEO  OATA  TO  FEB  TATA  ARRAY 

-C- - 

C  CHECK  FOR  BAD  CYCLE 


£9  9999*999 99 99 9999 9*999999 9-9 999999999999-99 99999999 99999999 999999*9-- 

C  CHECK  THAT  THE  PRESSURE  IS  IN  RANGE  9 

-£9  9999  *999 99 99  999999999  999*99999*999999*9*99*9 9999 99999999 99999999- 

c 

- It— LVAim  U6T.F0  OCW-l AO  . 

*  VARI  D.LT.-S)  GO  TO  999 

99  • 


C99999 9999 9999 9999 9*999*99 99 9999*9 999999 99 99999999 99999999*9999999 

C  - CHECK  THAT  THE  GATE--IS  IN  RANGE- — _ 

C8 999*  9 99999  *9 9999  99999 9999999 99999999 9999 999*9999 9*9999*9 99999999 

-C - 

IF  (VAR (INI .Gt • T POC W( 3t I .OR. 

- 8 - VAR(  lAl.LT  .FDOCW  (3  1)  )  -GO  -TO  999 - 


.,.6  ■■ 


F 

3  \ 


s 

L 


i 

i  >• 


II 

D 

I 

I 

6 
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C*  *****  *** •***  **  *• ********  **  *••***•*••«****•***•** ******** ******* . 

C  CHECK  THAT  THE  TI  ME  FIELDS  A  Ft  EACH  IN  Ra*GE  * 


. IF-  (VAR  (ISI  .CE.2A.QQ0  I  60  TO  9« - — 

IF  (VaR  Ui)  .GE.t-f^^CCn  I  CO  TO  90* 

_ _ If  (VAR (17) .GE.6r.3G0)  SQ_tO  9«* - 

IF  (!FIKST.hE.1>  CO  TO  *5 

C. _ 

C  COMPETE  PFT1"C  f°R  F  TRST  CYCLE  or  FIRST  S£6*CMT  ONLY. 

CPTIhE=VAr(  1H*v  Ah(i5  1/2*. * VAR(  1  *)/ 1**0.  «VAR(  17 1/86 ACC. 
JFlRSTrO 


-C — rCoHPWTt — R  EL-  StC— FO  R— ALL — C  yXL-E-^. - 

FELDAY=VAr( 1A)*vAr( 15)/2*.*VAR(lfcl/ 1990 . *V*R( 17 ) / 86  ACC. 


M2  =  *»2  ♦  I 

C - ORTHOGOM  LI-2E— V-ELOG-I-T-f— VtCT-OR - 

C 

- PO-99  10J=l-,3 - 

X0-0PTH06(1 ,j> 

nA  —  «  -v 


X0=V*R(  *ORTHOG|  I*  l,  J»*XO 


DO  60  1=1.3 

VMH4  *H)  =W  ( II - - - — — 

LENGTH  =  SQRT<w  ( 1  )*W(  n«W<?)  *U(  2>«V(  3)«U(  3!  ) 


FORMAT!  •  LENGTH  OF  ABSOLUTE  VELOCITY  PRIOR  TO  TRANSFORM 

— ,G20  ,S-> - 

WRITE  (8  ,2999)  U(  D.WI2I  *V(3I 

- FORM  Al(  i— V — — PEF  OR  E— U-r-3C20^7) - 

TO  150  1=1,13 

T*-t  _ _ _ 


V«( I,N2 >=V*R( I) 

-VN(-m,H21=0f>tlME - 

VH(15tN2)=RELSEC 


IF  (N2.LT.500.AND.NFLG.NE.1 )  GO  TO  10 


_  c-  — ***  *  *— - - - 

C  Write  FEB  SC6HENT 

C - - - — 

*0  18=0 

- if— (NfL6.CQ-.U-NW=N2 - 

*2*0 

_ 1100 CM  t«t  |=lbOCU(«)*l - - - — 

IOOCWU  3 ) -NBA  D 

- FOOcH(1C)=OPTImE - 

FOOCW(12)=FDOCV(  ?(.) >F DOCV (37) /2R.*FOOCV(  3*1/1990. 

_ S _ *F-aoC*.U9)/BKATn.  - - —  — 

CALL  ZVRTT  (2o,lF,IB| 

- 1  SEC  =  IS EG  *-l- - 

IF1NFLG  .£Q.  1)  RETURN 

. . KSAGza - 

C 


161 

C 

fO  200  7  =  1 .  IS 

• 

1  wc 

163 

ro  200  J=1 , 300 

27 

• 

1 


1 66  -  200  ■  —  TfWII.JIzO.OO - 

165  C 

166  - 999 - CONTINUE - - 

167  c  count  a t,n  omt  pao  cycles. 

168  r.BAD=N8  AO*l  — . . . . . 

169  NgOlOl  =  N8DT0T  ♦  1 

170  C  - - - 

171  - £*♦*♦*******♦****♦**•***************************•***•*♦***♦****•* ' 

172  C  WRITf  I Np OR  HAT i ON  CONCERNING  THE  BAD  RECORO  TO  UNIT  9  * 

l7]  C**************************************************************** 

176  C 

— 1-7-S - C - - - ^ - 

176  WRITE (9,80001  N2.ISEG .NBOTOT.V AR( II .VaR (161 ,V AP Cl  51 ,VAR( 16) 

177  8000- - F ORhA T I - 3 15 ,66 16,54 - - 

178  C  ***** 

179  C — TES  T  FOR-XNa-OF- -CAST - 

160  C 

- 181 - 1X-UOXO— -1— yANQ.— -K-F(,6— .f  0,-Q)  GOTO  IQ - 

162  IFINFlO.EQ.1  .AND.  KTLG  .EQ.  1>  G0T060 

185  - JF  -INFLG.NE  .1  )  GO  TO- 10 - 

186  C  IF  CEO  OF  CAST,  RETURN  TO  CALLING  ROUTINE 

185  C — AND  -TER Hi NA IE -FEB  XI LE. - - 

186  C 

-  187 - C — ***** - 

188  RETURN 

189  1020  WRITE  (6,1321)  L - - - - 

190  1021  FORMAT  (•  L=* .13,*  IN  DATABLOCX  READ.*! 

191  - END - - - 


FS3  =  rLo(Ct5tirSll)> 

■IF4F  $2- ,  NE  ,  -  2  5Sj  GOTO  i$ - 

1FIFS3.NC.31)  WRITE<16,87?7)  IFS(1)»FS3 


IF (F  S3  .NE .  31)  SOTO  1 5 
C  -***«* - - - 

c 

c  ***•*. - 

C  MOVE  NEXT  37  ByTCS  INTO  RO  ARRAY 

- MUttt  -  NH  B  Y-T£ - 

t-0  20  1=1,3  7 

- CALL  MO y£-< 1N3UFtNByTE,R0|X ),!,*> 

r;BYTr=NETTE4l 

20  _ CONTINUE - 


IF  (NM9yTE.CE.53A)  KFLS  =  1 
lFtNMBYTE.GE.S31!  )  NMBYTE 


CALL  DECODE  | RD  rvAR  > 
***** 


SPRTtS  V. DECODE 


iA*L 1BI1 ) .DCCoDE 


SUBROUTINE  DECODE  IRD.VAR) 


C  THIS  PROGRAM  decodes  the  data  broken  out 

#*  DV  CI'fiOAilTr  >.r  r\ 


C  ARE  STORED  IN  ARRAY  VaR. 


PEAL  R0(U,VAR(1) 

- IHTIGER-ItJ-U^) - 


L'«W 


COSI-2t»t2 - 

TII  )=FLD(0,8,r0I 1-1  )) 

ELo (2D,  A.T4  1 1 .)  =F  LOCO****© 1 144 


C  PRESSOR  E - 

C 

C****** *******  ********* *•***»< 
C***** ******** **•**•***••*•• *4 


********************* *********** 


PRESSURE  CALIBRATION  OF  IQ/8/S1  ,  10/9/81  APPLIED 


C**»«*  A***  ************************************************ *********** 


VAR(1)=  0.9989789*VAR(1 )*C. 20325 
-IF — (F LD-17,1  ,RD  (  9  )  1  ,£0.4  |  -VA* j+4*~VAR (  J.) 


C — SLOO-TC  KP. - 

VAR(2)=T(A)/20r0. 


conductivity  * - 

VAR(3l=T(6i/lC00. 

FAST  TCNR. 

- VARJ 4 1=4161/? COO- - 

IF  (FLO  (Stl  ,**0(9  II.EQ  .1 1  VAR(4I=-yAKl4l 

CONSTRUCT  VECTOR  COMPONENT  WORDS 


CO  IP  1=10,26,2 

- 1(  I -9l=.EL-l>T-C,  6,  R  r  ( Ii4 - 

FLO  122,  6, T I  l-9n=FLP|o,fl,R0(I*UI 
IF  (rU>(7,l,R0(lM.EQ.l)  HI-9lr-T(I-9» 

COMPUTE  VECTOR  COMPONENTS  -  - 


-VELOCITY  VECTOR - 


C***  - - - - -  *** 

~C*** — NOTE-  t-MO  CAirl-flR  A  7-lONS— HAVE— feEEN— ApR-tTEO—To-WE-  VELOCI-TV — *** 
C***  COMPONENTS  IN  INSTRUMENT  COORDINATES.  *** 

C*** - *** 


-VAIM5»=+f«  ♦rOO  61039 — 
VAR  (6  1=-T C3 I*. 0061039 
4>AR-t»=T<S)*.QG61039— 


VAR(M=-T(7l/l00C. 

- -VAR  (V  1— T( 91  /10CO  . - 

VaR(1C1=T(111/1C00. 


acceleration  vectop 

—  uinJ  «_i  1  11  It  /i  n  rn  - 

VaR(12»=-T(  1SI/1CC3 
- V*TM  1  *»-=£(  4-7 1/4  C  GO. 


VAR  (ml  =100 *FLO( 0,4, ROT  3711 «10*FL0< 4,4, RD(37)  ) 
- *FX0ia-r4-rR0(-36ll - 


C  HOUR _ 

VAR (1  SI =10*FLD(4,4,RO(36))*FLDfO*4«RC(3S)) 

r 


C  MINUTE 

- VAR-I16)-=10*FU)4  4,»rWM3S)  )«FL0(0,4rRD(34l  1— 

C 

-C — SECOND - 

VAR (1 71 =1C.*FL0 ( 4,4 ,RD( 34 ) I ♦FLO 1 0 ,4 ,®0f  3311 


♦FLDI4,4 ,RD(32))/10C0* 


RETURN 
END - 


30 


jA*LIBCl>.MOV£ - i. -  ■  -  .....  ■ 

|  SUBROUTINE  HOVE  cFROM tIF6VTC,TOvITs VTc.NRVTeSI 

- 2- - G - ; - 

S  C  THIS  ROUTINE  MOVES  BYTES  FROH'AM  INPUT  ARRAY 

- 4 - C-  -fBO*-,  TO-  .AS-QUTpU?  ARRAY- TO*..  If  BYTE  IS - 

$  C  THE  BYTE  NUMBER  OF  THg  fjRST  BYTE  TO  BE  TRANSFERRED 

- 4 - C — Fro  H— THE — INPUT— ARRAY-, — ITBylC — LG— THE- By  TC— NUMBER - 

7  C  OF  HE  POSITION  IN  THE  OUTPUT  ARRAY  WHICH  WILL 

- 8 - C — R£CEJVE-THE-TRAWSFER*ED-»YTE» — NryTeS-XS-ThE - 

9  C  NUMBER  o F  BYTEs  TO  BE  TRANSFERRED. 

- 10 - C - 

II  DIMENSION  FR9N(l|(T0m,MASK|l0ltl3H(9) 

—  ■  12 -  - C-AT  A — 18  v>f  28  .20.  1  2 .4  .0 . 2 4^  16  .8 . OF-  — — — 1 - —  ■ 

is  Data  mask/268435455  ,-<.67386680. -i04  4  4ao,-4oso, -is, 

- 14 - 4— 4278190C80  .-147116  8(3  , -6S2*C, - 

15  t  -255. 92  96 96  7295/ 

- 16 - CO-  1000-NHs  1.N8.YTES  . —  ■  ■  -  - 

IT  C 

18  - C.  .  •  «-*R£-T — A— BY-TE - 

19  C 

-  20 - NBSCNM-IHIFBYTE - 


21  INS CNB-lJ *8/36*1 

_ : — 22 - ICA4HBUSW0 *9 »C0.< NB-Sl  >■  60-  TO  100 - 

23  lBS-t  nB  -1  I *8-1 IN -1)  *36 

—  24 - IE — llBS^r  0.0)IB-S^Q - 

25  1NO=FlO(IBS,8,FROM{IN|) 

26  _ -60- JO  2  CO - 

27  100  IN0=FL0|32t4tFR0MliN)  l*i6*FL0C0,4»FR0HIIN*in 

- 2« - 200 - CqhIIwUE - 

29  C 

3q  , - C.  ..  ».PUT— A— By  TC. - 

31  C 

32  - _N8SIJIM^1)*I4BYTE - 

33  I8S=M00lNa-l,9l*l 

- 34 - WE  LSH  OP  ma-J^SLMm - 

35  IMS  1NB- 11*8/36*1 

36  - IF- Cl  NS  -S  7/9*9.  EQ„(  NS  *54  4-GA-TO-  303- - 

37  IN0=TN0*?**I3MIIPS» 

38  - T0(iN)=aN0IT0(In1.hASk(NELJ) - 

39  TOClNirORlTOClNl.INoi 

- 4Q - CO.  IQ  4 CO - 

41  300  T01INl=A:iDlT0llN»,MASMSl» 

42  - TO  1  IN  jr  OR  IT  01  IN  l  ,FL 01  28  ,4  » 160)1 - 

43  T0<IN*1  |  =  AN0(T0|IN*1).H«SKC10I) 

44  _ TO  ClN*H=0RCT0.1lN*l»-,FLDC32,4,IND»*?**32.1 - 

45  400  CONTINUE 

_ 46 _ C - 

47  1000  CONTINUE 

48  . RETURN _ 

49  END 


*rRT,S  V.MAPIth _ 

8A*L4B  1*4  *COM  VERSION - A - 

1  C  MAIN  ROUTINE  FOR  CONVERSION  PROGRAM. 

- 2- - e — : - : - 

3  c  THIS  program  converts  the  VELOCITY. 

- 4 - 6 AGCHE*-AT-10Ht»NP- MAGNETIC  VECTORS - 

5  C  FROM  INSTRUMENT  TO  6E0H A6NETIC  COORDINATES. 

— -6 - € - 


7 

- 8 


C  ******* ******** ******** ************************ ******* 


C  COMMONS  FOR  7 RE  AO 


COMMON/ ROATA/VRt ISf 5QCI 

-com  mon/roocf/fdo  co«4io-^— 

COMMON/ ROOC 1/ ID OCR 1 20 » 
-CONMOM/ROOC  A/  AOO  CR 1-504— 


- C - COMMCNS-FOR -CON  VERSION - 

COM MON/ BET A /DM A T (3, 3) .0RTH06 I 4 , 3 » ,A CCAL { 3 » 

_  com  nctcu .Accn »  r,m  ,u<3>  ,a  »3  ♦  M‘** 

- C — INI TIALI7E- VELOCITY  ORTMOGONALIz*TJON-M*TRIx - 

TATA  ORThoG(1,U  ,0RTH0GC7,l».0RTH0SlT,n  .0RTH05(4,l  j/ 

- *-_*.055,-,«E.J,l  .58  ,-.794/ -  - 

DATA  ORTHOC  Cl, 2) ,OSTHOG(2,2),ORThOg(3»7) ,0RTh03I4,2 1/ 

- r.0Ur.)5^ - - 

TATA  ORThOgI  1,3  >,0RTH0r,l2,3>.0RlH06l3.3>»0RTH0ol4,3»/ 
- 8_.6l2,.849,v.071,.oo»/ - - - 

C 

-  C-  INITJAL17E  ACCELEROMETER  CALIBRATION  COEFFICIENTS - - 

DATA  AC  CAL » 1> > ACCAL 12 », ACCALl 3)/-l. BOO, 1.030,1.000/ 


10  CONTTNUE 

- CALL -/READ -110, ir,0| 

C 


IPRC5  1Z*VL0CG1* 

- - 1PR  (6)= *VL0CG2< - 

JPR17  1='VL0CG3’ 

_ 1PR  (8  1=.*MaGB1  « - 

IPR  (*»=  *MAGG2» 

- IPS  II  Jl-tMABGI* - 

IPR (111 =* ACCLGl * 

_ IPR  112) =*  ACCLG2 * - 

IPR ( 1 31 -*  AC CLG3  * 

X — ADD  -I®  AOOc - 

AOOCP 13 7| :*  SEOMA» 

- : - AOOCR<3  8»-’-G-COOP’ - 

AOOCP  (39I-*  0  * 

_C - - 

ncycle=o 

c - 

20  CONTINUE 

- NCYCLE=NCYCLE>1 - 

IF  (MCYCLE.GT.MR1  60  TO  100 


50 - CONTINUE - 

C 

-  C  CAL  I _ *ECT  OR  CON VERSION -SUBROUT I NFS 


'x&miS&i--*- 


*  V  ‘  *'■  •  for a! 


32 


72 


CALL  VC  ICON 


-73- 

78 

75 

7b 

77 

78 
~J9- 


C  ***** 

C  -  WRl-7E-N£y-$EGMtN7- 
C 

too - contt  title 


IF  (K-CYCLE.LE.NR  )  GO  TO  20 
CALL  zw RJL3 -ImZO-t  IE-,61 - 


80 

81 

62 

83 

8*1 

8S- 


- C — C8E-CK— F-OR-  END-OF— CAST - 

IF  (IDOCRll l.NE. 1)  CO  TO  10 

- C - 

c  ******* ******************** ************ ******** ******* 

- E*D - 


8PRT.S  V •£ ONVV tCTO R 


i A*L  IB  (  l4-*CO»iV¥tCTOR - > - 

1  SUBROUTINE  VCTCON 

--  -2 - C - 

3  C  THIS  SUBROUTINE  CONVERTS  THE  VELOCITY.  ACCELERATION , 

-C - Atto  K»  6  Kr-I  I C-V  ext  pRS -F-gOW- -I  N  5  T  R  U«EW7-  -  CO  OBO  I  (HT&t - 


5 

b 

7 

9 

9 

-lO- 

11 

12 
13 
In- 
15 

-lb- 

17 

18 
1* 
20 
21 

-22- 


C 

-C 

C 

C 


INTO  GEOHASNETIC  cOOPoINATeS. 


************************************************** 

-COMMONS-  FOR— CON  VERSION- 


COMMON/ BETA /PmA  T  ( 3,  3»  .ORTHOG  l 9 , 3  1  ,A  CCAL  (  J 1 
-CON HON/ BOAT  A/  VB  MbrSO&l - 


C0HHON/RD0CF/FnOcR( 40 1 

-C0HH0N/R00C1/ID0CR(  20-1- 


— COMMON— NCVC  Lt,  ACCf34rXM  AOH->tUF31  «A  <  3)  ,M  «  3^ 
REAL  U  (3  I. LENGTH 


***************************************************** 


-*****- 


C 

-C- 


CALL  SUBROUTINE  TO  CREATE  TRANSFORM-MATRIX 


23 

2% 

25 

2b 

27 

28 


CALL  TRnnTX 


***** 


00 


COORDINA  TE 
-  TO  30  1-1- 


CONVERSIONS 
7 - 


29 

30 

31 
3? 
33 

Sk¬ 


ill  1  =  0 
-All  1=0- 


H(I  )r 0 

-00-30 -J=  1*3- 


3$ 

3b 

37 

38 

39 
-80- 


30 


U(I  )rU(  1|  «B  H»T(  I  .J)*VR(  j*i|,  NcYCLE  > 
illicit  I  LUBHATl  L^a*ltR<.Jt4S-.WCrCL£l- 


2  D  Cl 


1001- 


HI1)=H(  11*8  MATH  ,J1*VR( J*7,NC\CLEl 

- NRlTE(C*2001  I  Ul  lJ-.U(2  1.U131 - 

FORNATl*  V  -  AFTER  s’,3020.71 

- LENGTH=SQRT(UUl*U<l)*Ut2>*U<?>«Ut3>*Ul3)3 

write (9tioou  Length 
- EQBMA-TI*  ABSOLUTE-  VELOCITY  AFTER 


TEAL  SC  OiB-al  .G20.G//1 


81 

82 

83 

88 


-X — RUT.  NEW- VALUES— INTO  -OAT  A— ARRAY- 
00  SO  I =1,3 

- VR  (  I*  4.NC  YC  Le  l=U  (XI - 


33 


»***••**> 

RETURN 
END - 


*****••**- 


SPRTyS  V. TRANSFORMS 
3k*L IB<1 ». TRANSFORM 


SUBROUTINE  TRNMTx 


THIS  SUBROUTINE  CALITRATES  THE  ACCELERATION  VECTOR, 


AND  CREATES  THE  TRANSFORM  MATRIX  ‘eMAT*. 


******* ****************  ********************* ******** 

COMMONS  FOR  CONVERSION 

- - — COMMON/  BOAT  A/  MR  >  - 

CONMON/KOOCF/FDO  CM  I  AO  I 

- COHHON/fcE4A/PHA4  +  3,  3)-,0RTH06  *A^34-,ACCAL  (  2  I - 

- COMMON— NC-YCLE~» -AC  C ( 3  )  ,  X*A0+34^ U  ( 34-, A  +34  , H 1 3 1 - 


-C-  CAL  IFRA-TF— ACCEL-ERAT-XON-VCCTOR - 

PO  10  1=1,3 

-40 - VR  W*  10  ,NC-y  CLE4-C  ACCAL  ( I-f*  VR+ 1*1-0 rNC-VC  LE->  — 

C 


0 R HALIZE  *  A  - — AND  «  H  l  y EC40 


00  20  1=1,3 

xna  =jc  na_*vr  c  mo,  nct-cce-I-*-* 

XNH=XNH ♦VRl 1*7 ,N CYCLE )**2 


XNA=S0RTIXNA) 
XNHiSCR  T(*NHl 


ACC  (I!  =  VR<I«10,NCyCLC )/XNA 
-30 - xHAfcl 1)  sVR*<  I*  7,NCYCLt  J/XNH— 

c 

CALCULATE  D  IP_F  RQM_*  A»_  ANQ_»Hi — 
SlNDIP=a 

AO  SlNOIP=SINoIP»ACC<I »*XMaGCI» 

_ OlP=AslNtSllJ3IPl _ 

C0SDIP=C0SC DIP  1 


C  ***** 


c  crfate  transform  matrix 

_  00  50  1=1,3 

50  BMAI(3«1I=-ACC(II 


53  C 

54  .  _ XT2  =  l.0-EMAH2,l>**2-BNATC3,l)**2 - 

55  lr  IXT?.tT.oJ  *T?=-XT2 

56  _ PMUll.l)  =  5Qf»UJlT2» _ 

—57 - IF— 1XMA6T2I  .LT.G)  BHAHtrl >«-GHAT‘H-, H - 

Si  c 

—  59- - : - XT— ACCI-?l-*xHA6|-3  )“ACCl3)*XHAGf£) - 

60  BMATUt2t=-BHATll,l  >*UCC(1  >*XmA6(3  l-ACC  I  S»  *XHAS  1 1)»/XT 

61  - 1MAT<  1.  3>SpMxm  ,1)*<  ACCm*XHA6<8l"AeCt2>»*WA&  W-A MMf - 

62  C 

63  - C — ******************************************************* - 

64  C 

65  - RE-TURN - ' - 

66  ENO 


8PRT.S  V.MaPCON 


i  A*L  IB  4-1-1  -.PREFI X - > - 

2— - C -PROGRAM - v-,  PREFIX - 

3  C  PURPOSE  i  TO  1.  cOHPUTe  WI  FROM  PRESSURE  DERIVATIVES  AND 

- 4 - C - 2-..  ATTEMPT  TO  -EXX  ■Ef-FECTS  OF  0¥e*RA  ■  *M® - 1 

5  C 

6  - C**»*»*»*»6»** ***»♦♦♦♦♦»*♦*»♦♦♦»♦»♦♦*»***♦»*»»»»»» »**»»♦»»» 

7  C 

8  - C - 

9  C 

- AO - C - 


11  '  COMMON  /K9UFF /  K SEG, KNR, KREaL 

12  - COMMON/ BETA /UNORTHC  4,  3»-,ORTHOG(4,J»  ,wC3>  ,Q4-3) -  - 

13  C  COMMONS  FOR  ZRE AD 

14  - C0MM0N/RH0R/LR,NR,N6RtNNBR*NMFRfNFR-tNIRTNARtIPRll51  - — 

15  COHMON/ROAT A/VRt 15»500» 

- - 16 - COMNONAROOCF/F‘-&OC«140-» - - 

17  CONM0N/R0oCI/IDOCR(20> 

16  - COMMON/ RDOC  A/ ADOCRf  501 - - -  - 

19  C  ■  . 

20  - COMMON/DIA&$/N5GRtMSGNrNMNRfMNNU,NNAP,NNF,NNIrVNA,lRST,IHST 

21  C 

- 22 - C - COM  MOMS  FOR  ZURII - 

23  C 

24  - COMMON/ WHO(ULM-rNU,NBUfNMB»LrNMFN-rNF«,NI«fNAt|»IPut  151 - 

25  COMMON/UDaTa/VU( 15,5001 

26  - COMMON/ WD0CF/F00C«I4Q1 - - 

27  COHMON/ROOC I/IOOCWt 20  I 

C0HJKflN/MDQCA/A0QU(5flj _ 


29  C 

30 — - - C - - 

31  REAL  VROTt 5000) 

.  32 _ C - - - 

33  C***** ********  **************  **************************** ********* 

— 34 - c — initialize  .commons  fcr-zreao,  anil/wrjlu _ *. 

35  C***** ******** ************** ***•**•*’  ********************* ****••» 

36 

37 
3« 

3’ 

— 4a 

41 

42 

43 

44 


OAT  A  LR,NR,hFR,NIR,NAR/ 15 ,500,40,20  ,50/ 

- DATA.  L«  ,NN,NFH,NlW,NAU/AS,5C0,4Ci,23,50/  _  _ 

DATA  NNNR,NNN«,NNjP,NNFrNNl,NNA/503,500,15,43,20,50/ 

_ C - - - 

c***«* ************ ******** ******** ****** ******* ***********  ******* 

_  C  INITIALIZE  MATRIX  FOR  -CONVERSION-  BACK  TO  UNORTHOGONAL  ACOUSTIC* 
C  AXES  COORDINATES.  * 

_ c********* «*•*  **«*  *****  *********** •*******••*••••**«•••••«••*••*• 


35 


_ c _ 


r  at  a  lw  on  mil  »i  t-«mtoajtuz.».u->uKaaTH.i3.»4.>  ,uNoinmji  ,u/_ 

f  .534,  .02485  ,-.67493  6576q/ 

..  ..DATA  IW09Th(1,21  ,UN0RTHt2,2),UN0RTHt3,21,UN0RTH(4,2l/ 

S  .3505, .66162, .01352  , .63666/ 

.  CAT*  UUORTHU,lltUNORTH(2.3»,UNORTHl3t-31,UKORTHC4, 31/ 
t  .2268, .U3l98,.70c48 ,.60135/ 


C  INITIALISE  THE  ORTh^GOMLI  2*  TION  MATRIX 

C*M**M«M***,M«M*Ms'M*M»****M**»*«»M**,**M*,M*M*,* 

C _ - 

- EATA-ORTHOG  U  ,1  T-tORI^Os  4  l-lrOR7+K)fr<-3yl  l-yOR4HOS<-4-rl -7A 

*  -.055, -.803, 1.58, -.794/ 

- 0ATA-0RTH06  <1  ,2  »  ,0R  TH OG  <2 , 21  ,ORTHOG  (  3  ,2  »  <  THOS 1 4 , 2-»/~ 

*  -. 197, -.684 ,-.01 1,. 759/ 


*  .612, .849, -.371, .809/ 

- MSg  R: 1 - 

MSG  kr i 

C**MA*,**»********t6**m****,*M„***M*M,**,M#m****l*«M4**- 

C*  ENTER  INPUT  and  OUTPUT  FEB  FILE  INFORMATION  * 


MFITE  <  6,2001 - - 

200  FCRMaTI/*  ENTER:  NUlNl  ,NSCC  1  ,NSSEG  1*  1 

- - READ15,3001  -LUIN  1-,NSEG )  ,  NS SE  Gj - 

300  FORMAT  I) 

- Wf I T  C  ( 6,  4n01-WUT-N-HN6E<7l  .^5S-E-&i - 

400  F  CRM  ATI*  NUlHl  =  *,I2,  *  NSE61  =  • , 15 , •  NSSEFl 

- 700-  FORMaT|«-NUiM  =-*,12, NStG^-s-^-IS  , «-*SSE&2- 

*  •  MSGL12  =  * ,121 

—  C- - 

C  RE AO  SPECIFICATIONS  FOR  OUTPUT  FILE 


*,I5  1 

M5r 


WRIT  E ( 6, 9001 

-908 — FCRM*T<-«-  ENTER  NUOUT^-l - 

READ (5  ,3  00 1  NUOUT 

- WRIT  E ( 6, 9999 1  NUOUT- - 

9999  FORMAT (•  NUOUT  =  *,I2l 

_C - 

c  ******  *******  ******************  ***•*****«,*«•*  ********  ********* 

C*  BEGIN-  MAIN-SECTION-  OF  -PROGR  AM-HERE - *- 

C*  * 

C* - MI  MlLL-RtPLACE  THE  INITIAL-TIME  1VA»  -*l*  1 - * 

C*****  ********  ****************  ****************************  ***** 


NN  =  8 

_NX_i-2*NK,l  _ 

CELT z  0.0625 

K  =  NN  - - 

M  =  K*|K«1)*(2*K*1) 


a  half  midTh  of  derivative  smoc 

-8  SET  SIZE  OF  data  wInOow - -i 

8  NOMINAL  TIME  D1FF EPENCE  . 


8  Op/DT  T.  HI  CALCULATION. 

CAU.  2BUEFU1  ,S3,NUIN1,NX,nSSEG1  ,VROTl - 

IFIICOCR (l i.EO.l lIDCCwil 1=1  8  $El  ENO  Of  CAST  FLAB  IF  RE) . 


-T>Z_OU TOUT  0 CC-. — BLOC 


00  70 CC  IPrl  ,  40 

FO  OCk  1  IP  )=F-OOCR<  IP  1  - 


109 

109 

110 
111 
11? 
113 
119 

115 

116 

117 

118 

119 

120 
121 
122 
123 

-12ft- 


7000  CONTTMIE 
_ DO— 13C1—1P-1  «2£- 


7001 

7002 
-7qo3- 


ILOCu  IIP) =IDOCF|IP> 
00  7002  IP=  1,60  - 

AOOCW  |IP)=A00CPIIP) 
00  7003  IP=1,1S  -  - 

_I  PW4I  P4-^I  PR  1 1  P-1 - 

IPWIIK)  =  »«I  • 

-N  )i*  W  -=-NM  gR - 

NhFu  =  NHFR 


NPOINT  r  NSe0l*S00 


•MAXIMUM  POSS.  NO.  OF  OAT*  CYCLES 


C - 

C 

C*  ♦♦***♦**  ********  *****-*******************  •••***»•  ********  ******* 
C***** ******** ********* ******* ************ ******** ********  ******* 
-Xftftft _  *** 


125 

126 

127 

128 
129 

-130- 

131 

132 

133 
13ft 
135 
436_ 
137 
138- 
139 
lftO 
1*11 

-192- 
19  3 
199 

195 

196 

197 
-19ft- 

199 

150 

151 

152 

153 
159 
155 
ISO 

157 

158 

159 
160- 
161 
162 
163 
169 

165 

166 

167 

168 

169 

170 


C*** 

C**« 


BEGIN  THE  MAIN  LOOPS 


*** 

-***■ 


c* ****************************************************  *********** 

-  C* ********************************************************  ***«♦**- 

c 

_x _ 

00  1000  I  =  1  ,  NSEGl 

- 00-3000  1 1 =-1,5 CO - - - 

00  3000  JJ  =  1,15 

-VW|JJ,II|  -=  WR-I JJ, II  ) - 


3000 


CONTINUE 
- 00-1001- 


-j-j-i,  son 


ICVC  =  II-1)95C0»J 
-X-£4— I  CYC  -.LE.— NN*  l )- 4FTR=-ICyC- 


I FI  ICyC  .GT.  NNM)  IPTRr  NN* 1 

C***** ****************************************************  ******** 

_ C - C-CM.PU.TE-  M.I—AT  .  IHf  I-tJU-POINT  l-APPgQXUmH-Y->  _P LA  CE—LN - *_ 

• 

:  *6 »**» *****  ******* • 


c  variable  19  in  the  output  buffer 

c* ****************  *************************** 

c 

-OPOT  =  0- 


00  2000  K  =  1 ,NN 
- L£=MN.*4.*K _ 


2000 


LM  =NN*1-K 

-  o£DT=DPDT*K* IVRXI 1,LP>-VRX (1 ,LM) ) 

CONTINUE 

_OPOT=FACTOR*DP0T  - 

V HI  19  , J )  =  OPOT 


C****< 


C  —  TEST  FOR  GROSS  INEQUALITY  OF  Ul  AND  V.  ATTEMPT  TO  -  • 

C  CORRECT  *  BY  * 

C - 1.  RACK  TRANSFORMING  U,V,W  To  01  »Q 2 1 0 3 - • 

C  2,  REPLACING  Pj  ANO  03  BY  .657WI-.539  AND  * 

-C - .60185V1-.227,  RESPECTIVELY— -AND - *- 


3.  Rc  TRANSFORM  IN $  THC  Q  COMPONENTS  TO  U»V,*NO  b* 


IFI  A BS IDPDT-V RXI 7, 1PTR ) | .LT.42 ) 

1  GOTO  5000  8  TFST  FOR  PAD  VEPT.VFL. 

c* ****************************************  ***************** - 

c*  convert  velocities  back  to  original  acoustic  cooros  * 


-Dft  2500  K=l,3- 


37 


V  =  *♦«< JJ»*UNORTH( 

- CONTl NUe - 

fi(K)  r  V 

- CONTINUE - 

-Q4-H— s  U»6St»0RUT*-  0.»» -4  CORRECT  I-ST-ACOUSTTG-ATUS- 
0  431  =  0«  602*0P0T*  0.227  8  CORRECT  3-RD  ACOUSTIC  AXIS 


C - RE  TRANSFORM— RACK-70— ORTHO CON AL — COORD INA  TES - - 

C****« ********************** **** ********************** *»****•»■ 

-C - 

00  2610  KK-  1,3 

- U(KK)  -i -ORTH  06  (1  ,K-K  ) - 

00  2611  LL=  1,3 

- W<KK-»_=--y<KK7-*_0<LL  »*0RTH05«LL«1,ICK4 - 

2611  CONTINUE 


2610 
.5000 - 


V  V  <  KK  «q  ,J 
CONTINUE 
-Conti  nuE-- 


lt  XM- IF-U  -IS -CO on 


1000 - CONTINUE- 

STOP 

1 - CONTINUE— 

C 


IF4ICYC  .GT.  nn  .OR.  I  .61.  1>  C»LL  S TEpX If  1 , %2 , 1 1 
- CON  TINUE. - aiCOCI  1 1=1  ON  RE  A3  UNI  T _ 

continue  send  of  segment  locp 

-I-FII.  EC^uEXGi-1100 CU  Ll-lJ-i-  a- $E4-  EmD-CAST -F LAC--CN  -  L AS  T— Ss 
CALL  Z«RITInUoUT,IF,C  »  Si  WRITE  SEGMENT  TO  OUTPUT  FILE 
MTINUE - -B-  CND  0F-NA4N-  PROCESSING  LOOP 


CALL  IWRI TlNUOUTt IF ,01 

_ STOP _ _ _ 

END 


»EI»F  RETURN  On  READ  UNIT 

— B— SE4— E-MO-OR-  C AST  -FLAG - 

»  WRITE  LAST  BLOCK 


aPRT,S-V.2BUFC - 

ODE’T  1*FCHFILE  1 1 1  )  .VFIX1 


THIS  PROGRAM  COMPUTES  THE  PROFILER  CORRECTED  VERTICAL 

velocity  using  integrated  accelerations 
before  executing,  the  input  and  output  feb  files 

MUST  BE  ASSX6nED  TO  SOME  UNIT  NUMBERS 
MAPPING  ELEMENT  I  FCHFILEl .MVFIX1 

DIMENSION  A(1,5|,X|1  >,WI  (SC0),AI«I1  IN  ,S)  tDUtSO0),DWTIHE(5OO) 

COMMON  /RHDR/LR,NR,NBR,NMBR,NMFRVNFR,NIR,NAR(IPR(20I 
COMMON  /ROaTa/VP(IOOOO| 

COMMON  /RDOCF/FOOCRIROj 
COMMON  /ROOC1/1OOCRI20) 

COMMON  /ROOCA/aOOcRIZOOI 

COMMON  /DIAGS/MSGR,hSGw,NNNR,NNNW,NNlP,NNF,NNI,NN»,IRST,IWST 
COMMON  /OATA/A,Tl,NUMT,NUMP,NUMV,lSTART,TLAST,ALAST»PRESSltWl»FltX 
COMMON  /0ATA1/AIJ1,WSA VC,TSAVC,ASA VE,TlMDIFtINr.  X,GRAV,IFTRST,NUM* 
COMMON  /DATA2/WCHG,0W,DwTlME  ,I0w,F2,yPAR,NSSEG,.N0EXS,IUNR,IU'IWtM 
DATA  LR, NR, NpR,N i R,NAR  720,500,10,20, SO/ 

OATA  NNNR,NNN«,NNIP,NNF,NNI,NNA/500, 500, 20,10,20, 200/ 


c 


GRAV-O. 

A  CORK =9.  9« 

R HOG  1?. 9945 


[ 


I. 

i. 

i: 


i: 

i: 

i: 

r. 

E 

E 

[ 

I 

I 

I 


21 

22 

23 

24 

25 


29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


54 

55 

56 
«T 


60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
•1 
82 


C 

C 

c 

c 

c 

c 


c 

c 


2 

3 

C 


4 


C 

C 

c 


5 

C 

8 

C 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


6A 9=  MAXIMUM  TIME  (SECONDS)  BEFORE  SYSTEM  IS  SOLVED 

and  corrected  values  are  confuted 
SAFI =  MAXIMUM  TIME  Sa*  (SecONoS)  before  linear 

INTERPOLATION  FOR  HISSIN6  DATA  IS  PERFORMED 


6AP=1. 

6  API  =0.1 

IFIrsTzO 

WRITE  <  6.  1 )  --  - 

FORM  ATI*  INPUT  UniT  NUMBER,  NUMBER  OF  SEGMENTS.  STARTING  SEGMENT 
f HESS  AGE  LEVEL*. 

READ  <5  ,2  IIUNR,NSE6tNSSE6 ,HSGR 
FORMAT  (I 
WRITE  (6  ,3} 

format <*  output  unit  number,  message  level*) 

READ  <5,2  HUnW.mSGW 

NINO  Ex  =0 

IBGN =nsseg 

00  100  ISEG=1,NSE6 

call  2Reaoiiunr.if»nsscg) 

IFI1F.nE.0IG0  TO  99 
iflag=o 

NNSSE6=NSSEG 

LOOP  DETERMINES  VARIABLE  LOCATIONS 
00  5  J=1,LR 

IF  <IPR < J ) .E  0. *  VLOC 1 3  *  >  NUM V= J 
IF  (I PR <JI.E0.*ACCLI3*)NUMAZJ 
IF llPB IJI.EC.’PRESS* )NUnP=J 
IF IIPR (J| .EO.*  RELSEC* ) NUMTzJ 
continue 

s 

INDEX rNlNOEx 

istart=indcx*i  . . 

CALL  SUBROUTINE  TO  INITIALIZE  STARTIN6  VALUES  FOR  THE 
FIRST  REAO  OR  FOR  THE  FIRST  READ  AFTER  A  TIME  GAP 

IFIIFIRST.EO.OICALL  INIT 

00  10  I-IST  ART  »NR 
IN0=(I-1)*LR 

n  MOIF  =V  R I 1 No«NUMT I -TL AST 

IF  TIME  DIFFERENCE  TOO  LARGE.  SOLVE  TME  SYSTEM  * 

IFITIN0IF.6T.GAP 160  TO  11 

IN0£X=IN0EX«1 

ACCrVR  <INO*NUMA)  *ACORR 

IF  DATA  VALUES  ARE  NISSIN6.  PERFORM  INTERPOLATION 
IF  <TlM0IF.6E.6AP il CALL  CORREC< IND.MLAST.GAP) 


39 


E. 


s 


c 


•  9 

C 

IF  NO  DATA  VALUES  MISSIN6,  PERFORM  ACCELERATION 

86 

c 

integration  (Fit 

87 

c 

U 

!FITlN0IF.LT.6APi)FirFi*(VR|IN0+NUMA>*AC0RR*ALAST>*TlM0IF/2. 

•9 

c 

90 

c 

SAVE  LAST  VALUES  FOR  NEXT  CYCLE 

91 

c 

92 

ILAST=VR|INO*NUHTl 

93 

... 

ALAST=VR(IND«NUMA)«AC0RR 

99 

MLAST=VR(IND*NU«V» 

95 

PLAST=VR(INO«NUMP) 

96 

c 

97 

c 

T2=  TOTAL  TINE 

98 

c 

Hr  MEASURED  VELOCITY  (CONVERTED  TO  NETERS/SECi 

99 

c 

WfiA«  =  AVERAGE  VELOCITY  (DASEO  on  pressure  changei 

100 

c 

101 

F2=TL*  ST-Tl 

102 

U=GRAV*F2«Vr(InO«NUHV)/100 

103 

wbaR=o. 

109 

IF IF2. NE.O. IUB AR=- <PR ESSl-VR (IND«NUnP) )*rh0G1/F2 

105 

c 

106 

c 

SUN  NEW  VALUES  TO  ARRAys  IF  NOT  ALREADY  DONE  THROUGH 

107 

c 

The  correc  subroutine 

1«* 

c 

0* 

IFITINDIF.LT.GAPdCALL  ARRAYS 

in 

10 

continue 

112 

Go  TO  90 

113 

c 

I**." 

e 

r 

PR06RAH  CONTINUES  HERE  IF  A  TIME  GAP  Is  DISCOVERED 

»- 

V 

11 

NlNOEx=lNOEx  f 

117 

c 

118 

c 

TAKE  CARE  OF  THE  CASE  UHERE  A  TINE  GAP  OCCURS  BETWEEN  SEGMENTS 

119 

c 

120 

I F  <  IN  DC  X  .  EQ  .  0 1 NNS  SE  G7NN  SS  EG-I 

121 

IF* index. EQ.OJNINDE X=nR 

122 

c 

123 

c 

CALL  SUBROUTINES  TO  SOLVE  THE  SYSTEM  AND  COMPUTE  CORRECTIONS 

129 

c 

125 

CALL  SOLVE(PLAST) 

126 

CALL  WICOMPl AC ORR  » NINDEX# NNSSEG f IBGN t IFLAG, -API | 

127 

c 

12® 

90 

NSSE6sNSSEG*1 

129 

c 

130 

c 

IF  SYSTEM  HAS  NOT  BEEN  SOLVED,  RESET  INDEX  ANt *  SET  FLAG 

131 

c 

FOR  VARIABLE  INITIATION 

132 

c 

133 

IF(lFLA6.EQ.O)NINDExrO 

139 

IF<IFlA6.Eq.0>1FIRST=1 

135 

c 

136 

c 

IF  A  COMPLETE  SEGMENT  HAS  NOT  BEEN  REAO,  RETURN  TO  REAO 

137 

c 

WITHOUT  INCREMENTING  THE  LOOP  COUNTER 

138 

.....  c 

139 

IFt INDEx.LT ,NR }60  TO  9 

100 

CONTI  NUE 

192 

c 

193 

c 

SOLVE  SYSTEM  ANO  COMPUTE  CORRECTIONS  AFTER  LAST  READ 

199  ' 

.  c  - 

195 

NINO£X=INOEX 

196 

nnsseg=nsseg-i 

40 


14T  CALL  SOLVE(PLaST) 

148  CALL  NIC0MPIAC©RR,NIN0Ex,NNSSEG,IBGN,IFLAG,6aP1) 

149  $0  TO  999 

150  99  NRITEI6, 7) 

lsi  7  format </,•  read  error*  » 

152  999  E«0 

OOr*T19rCHFUElll  ).1MiT/¥FIX1 

C  THIS  SUBROUTINE  INITIATES  OaTa  FOR  The  BEGINNING  OF  THE 
C  PROGRAM  OR  AFTER  A  TINE  GAP  HAS  OCCURRED 

3  C 

9  SUBROUTINE  ZNIT 

5  DIMENSION  A(4tS>,Vlt500),AIJit4,S)vXt4>»DlHSOO)*DUT!NEt5tm> 

V  ‘  COMMON  /RDATA/VR I  100001 

7  COMMON  /RHDRsLR'NR,NBRtN*BRtNMFR,MFRtNl**NARtlPRt20! 

8  COMMON  /OATA/A.Tl ,NUMT  ,NUMP,NUMV,  1ST  ART.TLAST,  aLaST.PReSS  1 ,111,  FI «X 

9  COMMON  /OATA  1/ AlallyNSAVE  (TSAVE  (AsA yE  ,TlMOTF, INDEX  ,GR AV,IFIRST,NUMA 

10  COMMON  /DATA2/NCH6,DW,DUTXHe,IdW,f2,WBAR,N$SE6,IN0EXS,1UnR«IUnU,N 

11  C 

12  INDE-xSsINOEX+l  . - 

13  I0N=0 

14  IND=(ISTART-1)«LR 

15  C 

16  C  T 1(PRE SS 1=  START  TIME  AND  PRESSURE  FOR  INITIAL  READ  OR  AFTER  3A» 

IT  C 

18  Tl-VRtlNO+NUMTt 

19  TLAS  T=T1 

20  PRESS1=VR<IND*NUMP) 

21  F1=0 . 

22  MCHG=0. 

23  A LAST-0* 

24  VSAVE=0* 

25  TSAv£=T1 

'  '  ASAVE=0. 

C 

C  LOOP  ZEROES  ARRAY  VALUES 

29  C 

30  . DO  10  1=1,4 

31  00  10  «!=1  ,5 

32  AX *11  (I  ,J)rO  * 

33  10  A( 1, J> =0, 

34  return 

35  END 
«0Fm«FCHFILEl(l  »  .CORREC /VFI  xl 

c  This  subroutine  inperpolatEs  ano  fills -.gaps  for  missing  data 

c 

3  SUBROUTINE  CORREC  (IND,nLAST, GAP) 

4  COMMON  / RHOR /LRfNRf NBR ,NMBr, NHFR ,nFr ,NlR »NAr ,IPR (20) 

s  COMMON  /ROaTa/VRUOOOO) 

6  '  COMMON  /0ATA/A,T1 ,NUMT(NUMP,NUNV,ISTART,TLaST,aLaST,PRcSS1,MI,F1,X 

7  COMMON  /OA TA 1/ AIJl, HSAVE,TSAvE, ASA VE.TIMD IF, INDEX t6PAV,IFTRST»NUMA 

8  COMMON  /DATA2/MCHGfDN,DVTlME,XDW,F2,NPAR,NSSEB,IN0EXS,IUKR(IU9U,N 

9  DIMENSION  Al 4,51 ,NI« 5001 ,A 1.11(4 ,5 | ,X 1* ),ON(500 ),DVTI ME (SCO) 

10  C 

11  0U1=VR«IN0«NUHV1 -BLAST 

12  SLOPErOul/TlMOIF 

13  C 

14  C  COMPUTE  THE  MAXIMUM  NUMBER  OF  TIMES  TO  LOOP  IN  OROER  TO  DETERMINE  THI 

15  C  NUMBER  OF  DATA  POINTS  MISSING  IBASCO  ON  AN  OPTIMUM 

16  C  TIME  DIFFERENCE  OF  0.0625  SECONOS) 

IT  C 

18  LOOP=INT (GAP /O *0625 ) *1 


41 


»’ 

chkpts=o 

20 

C 

21 

C 

loop  determines  the  number  of  POINTS  MISSING 

22 

c 

23 

do  10  I=1,LOOP 

2% 

'  ■  "" 

PTS=AB$(TIHoiF /I -0*0*251 

25 

IPIPTS.LE .CHKPTS)  NP  TS=I 

10 

c 

CMXpT  S=PT  S 

*.6 

V* 

c 

01 VIDE  THE  TIME  INCREMENT  EOUALLY  A*0  COMPUTE  THE 

29 

c 

velocity  increase  per  increment 

30 

c 

31 

T INO  IF  sT  IM0I F/NP  T  S 

32 

OW2=SLOPE*TlMOIF 

33 

c 

35 

c 

LOOP  FINOS  INTERPOLATED  VALUES  OF  Fi ,  U,  E2»  AND  SUMS  TO  ARRAYS 

35 

c 

36 

.... 

DO  20  I=1,NPTS  . .  . . . 

37 

F l=Fl ♦DW2/1Q0 

38 

U=IHlAST«DH2«I  1/100 

3* 

F2=F2*TIMDIF 

50 

CALL  ARRAYS 

61 

20 

continue 

62 

C 

. .  . . 

63 

C 

RESET  TIME  DIFFERENCE 

66 

C 

65 

TIM01F=TIHOIF*NPTS 

66 

return 

67 

END 

OOFMl*FCHFILCl(ll.ARRAVS/VFIXl  . 

c 

THIS  SUBROUTINE  COMPUTES  VALUES  IN  THE  MATRIX  u$EO  TO 

c 

SOLVE  FOR  ALPHA,  BETA,  GAMMA,  AND  LAMBDA 

3 

c 

6 

subroutine  arrays 

5 

DIMENSION  A«6,SI,Vl(500l,AIJll6*5)»AIJ2(6,5l,DHIS001,DUTlME(530) 

6 

DIMENSION  X<6) 

7 

COMMON  /OATA/A.Tl ,nUmT ,NUMP ,NUNV, 1ST ART, TLAS T,ALAST, PRESS  1 «Hl, Fl ,X 

8 

COMMON  /D«T«1/aIJ1,USaVE,TSAYE,ASAVE,TIWDIF,InDEX,GRAV,IFIRST,NJMA 

9 

common  /data2/vChg,dn,ovtime,idw,F2,wbar,nsse6,inde*s*iunr,iunn,n 

10 

c 

11 

c 

ARRAY  AIU2  STORES  VALUES  FOR  THIS  DATA  POINT 

12 

c 

APRA.Y  AIJ1  STORES  VALUES  FOR  THE  LAST  DATA  POINT 

13 

c 

16 

• 

AI42U  ,1  l=Fl*Fi 

IS 

A  1*12  <1  ,2  l-Fl 

»’* 

A1J2U,3I=-F1*F2 

17 

A1J2 ( 1  ,6 l=Fl/2 

18 

AIJ2U,5»=Fl*N 

19 

AIJ2(2,1UF1 

20 

AIJ2  (2,21=1. 

21 

AIJ2 (2,3 1=“F2 

22 

AIJ2  (2,6  1=0*5 

23 

AIJ2  (2 ,5  )=N 

26 

AIJ2(3,1  »=Fl*F2 

25 

A1J2 (3,2 1=F2 

A1J2 (3 ,3 1=-F2*F2 

AIJ2 (3 ,6 )=F2/2 

A1J2 (3 ,5 l=F2 *V 

29 

AIJ2  (6  ,1 1=F1 

30 

AIJ2 (6  ,2 1=1* 

31 

A  1*12  16  ,3  1=“F2 

32 

AIJ2 (6 ,6  1=0* 

42 


*IJ2 (4 .5 )=MB AR.GR AV*F2 


35  c  loop  Integra  tes  ai<I2  values  into  array  a*  and  stores 

36  C  A1J2  VALUES  INTO  ARRAY  AI.'l 

37  c 

38  *  DC  10  1=1,4 

3V  DO  10  <1=1*5 

40  AII,JI=A(I,J|*(AIJIC I , JI+AI J2( I,  J) I •TINOIF/2 

41  AIJ1II ,J)=AIJ2(I,J> 

42  10  CONTINUE 

43  A|4*SI=W8AR*F2 

44  RETURN 

45  END 

FRT«S  FCHFILEl.SOLVE/vriXlt.WICOMP/Vf IXl,.WRlTE/VFIxl 
<9E331<FCHFILE1<1 )  .SOLVE/ VFIX 1 

C  THIS  SUBROUTINE  SOLVES  THE  MATRIX  EQUATION  AND  OUTPUTS 

c  start  times  and  end  times  for  time  and  pressure  for 

3  C  THE  PERIOD  SINcC  THe  LAST  TlM£  SAP*  ALSO  OUTPUT  ARE 

4  C  •  alpha,  beta,  GAMMA,  and  lambda 

SC 

SUBROUTINE  SOLVE  I PLA  ST  » 

7  DIMENSION  AA 14 ,5  I ,A*4 , SI ,X U » ,MI* 500 1 ,AIJl (4 ,51 

8  COMMON  /RH0R/LR,NR,NBR,NMBR(NHFR,NFR,NIR,NaR,IPR*20) 

9  COMMON  /ROATA/VRI 10C00 I 

10  COMMON  /DATA /A ,T1 ,NUMT (NUMp,NUMV, I ST ART ,TLAS  T, AL AST* PRESS  1 »UI*F1 *X 

11  Common  /datai/aiji,usave,tsave,asave(timdif,inoex,6Rav,1first,njma 

12  C 

13  NRITEI6,S1IT1,TLAST,PRESS1,PLAST 

14  51  FORMAT**  start  Ti»«£,  END  TIME  *,2F20.4,/, 

15  •  *  START  PRESSURE,  CND  PRESSURE*, 2F20. 4, /I 

16  oo  20  1=1,4 

17  00  20  <1=1  ,5 

18  20  AA (1 ,J Ir AC  I ,J) 

19  XMAT=SIMUl<4,AA,X,0,0,9I 

20  URITE(6,14)XMaT 

21  m  format*'  xkat  =',gjo.s» 

22  HR1TE*6,6)X( 1),X(2>,X<3>,XU) 

23  6  FORMAT**  ALPHA  =* ,610.5, 3X,*BETA  =*,610.5 ,3X , 'GAMMA  =*,Sxo.5, 

24  «3X, ’LAMBDA  =*,610*5) 

25  RETURN 

<  END 

0OF'3i*FCNFILEl(l  l.UICOHP/VFIXl 

C  THIS  SUBROUTINE  COMPUTES  ThE  CORRECTIONS  TO  VELOCITY 

C 

3  SUBROUTINE  wIC0HPtAC0RR,NlNDEX,NNSSC6,lBGN,.fLA6»GAP  l) 

4  DIMENSION  A*4tS)fX|4),MII500l,AIJlU,5>,0tf(500),0wTlNei5Q0) 

5  COMMON  /RHDR/LR,NR,KbR,NMbR,NHfR,NfR,NiR,NAR,IP' 

6  .  COMMON  /ROAT A/ VR 1 10000 ) 

7  COMMON  /DATA/A, T1,NUMT(NUMP,NUMV*IST ART, TLAST,Al  * ST,PRESS1 ,Wl,Fl  ,X 

•  common  /datai/aiji,nsave,tsavc,asave,timoif,index,g'  V,IFIRST,N'JmA 

*  Common  /0ATA2/WCH6fDy,DvTlH£,I0t|,F2,B8AR,NSSE6,lNDEXS,IUNR,IU9M,W 

10  C 

11  ISTARTsINOExS 

12  IFIRST=0 

13  IEN0=NNSS£6-ISGN«1 

14  DO  10  WRITE  =l,l£ND 

15  C 

16  C  REAO  SEGMENTS  BEGINNING  WITH  START  SEGMENT  OR  THE  .AST 

17  C  SEGMENT  WHERE  A  TIME  GAp  OCCURRED 


. 

t . 

4 

■■■■■ 

r 


L  f  i 


19  CALL  ZrEA01IuNR,IF,IBGNI 

20  C 

21  C  SET  LAST  CYCLE  To  RE AO  '  * 

22  C 

23  I  NO  Ex  I=NR-1 

29  1F(  III  RITE  .EQ.IEND ) INDEX l=NINDEx-l 

25  C 

INO-t ISTArT-1 ) *LR 
CH6TIH=VR  fINO*NUMT)  -TSAVE 

t*  C 

29  C  COMPUTE  INSTRUMENT  VELOCITY  tWXl  FOR  THE  FIRST  CYCLE 

30  C  IVX  VALUES  ARE  COMPUTEO  BY  EITHER  INTEGRATION  OF 

31  C  ACCELERATION  OR  BY  CHAN6E  IN  VELOCITY,  DEPENDING 

32  C  ON  THE  TIME  GAP) 

33  C 

39  IFf CNGTIM  »GE»6AP1 )VI( ISTART) =WS AVE«XI 1 J*( VR fIN04NUPV J -WLAST S> / 

35  8  1C0-X t3)*CH6TlM 

36  IFICHGTlM.LT.GAPllMHISTART»=WSAvE*IXtl)*IASAVE*VR«lND» 

37  •  NUMA)«2*AC08R)-2*<gRaV«X(3)>>*CHGTIM/2 

3B  ~  C 

39  C  FOR  FIRST  CYCLE,  WlzBETA 

90  C 

91  IFf IF IRST  *EQ«0 )Vl ( I  ST  ART ) =X ( 2) 

92  C 

93  C  COMPUTE  Ml  VALUES  TO  END  OF  SE6MENT  OR  TO  LAST  CYCLE  READ 

99  C 

95  00  30  Irl START ,INDEX1 

96  lNO=II-ll«LR 

97  CHGTIh  =VR  f I N0*LR  4NUMT I “VR  f lNO*NUMT ) 

98  IFfCHGTIM.GE.GAPllwHI4l|=WltH4Xfll*f  VRfIND4LR*NUMV)« 

99  •  VRflNQ*NUMVn/lt3-Xf3>*CH6TIN 

SO  IFtCHGTIH.LT.G«PllkT(I«ll=Vl(X  I «  «  X  (  11*  CVRI  lNo«NUMAI  ♦ 

M  >  VR<IND*LR*NUMAI«2*AC0RRI-2*(GRAV*X(3>) >*CHGTIM/2 

30  CONTINUE 

C 

59  INDEX 2“INDEXl*l 

55  C 

56  C  CALL  SUBROUTINE  TO  WRITE  CORRECTED  VELOCITY  TO  OUTPUT  FILE 

*7  C  "  '  .  ~  '  - 

CALL  wRlTEf IF , ULASTS, INDE X2) 

C 

60  C  RESET  VALUES  For  reaoing  next  segment 

61  c 

62  XB6NzIB6N41 

63  IFIRSTrl  ~  . .  . 

69  ISTART=1 

65  10  CONTINUE 

66  C 

67  C  SET  VALUES  FOR  RETURN  TO  MAIN  PR06RAH 

68  C 

69  IFflNDEX  *NE* 0)NSSEG=NN SSEG-1 

TO  IFflNDEX .NE.0)IB6N=IB6N“1 

71  IFLAGri 

72  IFIRST^D 

73  IFflNDEX .EQ.OlNlNOEXro 

79  RETURN 

75  ENO  "  ’  . . -  '  ■ 

«0FMI9FCHFlLElfll.WRlTE/vFIXl  - 

c  this  subroutine  writes  corrected  velocity  values  to  a  new  feb  file 

C  I  THE  output  FILE  CONTAINS  2  VARIABLES  more  Than  the  input  file) 

3  c 


I 


r 

r 

t 


4 

SUBROUTINE  WRITE < IF, WLASTS ,IK0EX2 1 

5 

CBmnOn  /RHOR/LR»NRtNBR «NMSR,NMFR,NFR ,NIR,NAR ,XPR (201 

6 

COMMON  /ROAlA/VRC 10000 > 

7 

C CNN ON  / ROOC  F/FOOCR 1 40  I 

• 

COMMON  / RDOC 1/ 100 CR ( 20 1 

9 

COMMON  /RD0Ca/ad0c4(200I 

10 

common  /whor/lw«nw»nbw ,nmbw«  nmfv*nfw »nin , nan ,ipw (201 

11 

COMMON  /W0ATA/VW1 10000 1 

12 

- 

C  OHMON^/nPOcF/fdO  CN( 40 1FW00C I/IOOCWC 201/NOOCA/AOOCWt 200* 

13 

COMMON  /0IAGs/*s6R,HS6w«NNNR»NNNtf»NNZP,NNF,NNXtNNA(XR3T,XWST 

19 

COMMON  /0ATA/A,Tl,NUMT,NUMPtNUM4.ISTART,TLAST,ALAST,PRESSl,Ml,Fi,X 

IS 

COMMON  /0A1Ai/AIJl,WSAVt , TS A VE t AS  A VE »TlMo IF, INOEX ,6ft  A V,IF iRSTt NJMA 

16 

COMMON  /0AfA2/NCH6,Dw*0wTXME,X0WtF2*wBAR,NSSe6,IN0EXS,XUNR,IU4w,H 

17 

DIMENSION  AU,S».WH500»  .A  I  JlU.5  »  .0  WC500  »  tOWTlME  I  500  ) 

18 

C 

1* 

00  10  I=1»NNF 

20 

10 

FDOCW  (X  )=FDOCR  (XI 

21 

00  20  1= 1 , NN 1 

22 

20 

XOOCU  ll  l>  IOOCft  (It 

23 

00  30  1=1, NNA 

24 

30 

AOOCW (l)=AD0CR(Xt 

25 

NW=NR 

NRUrNBR 

nmbw=nmbr 

BL0 

nmfh=nmfr 

24 

NFwzNFR 

30 

NIU=NIR 

31 

NAU=NAft 

32 

Lv=LR*2 

33 

NUMV  lrNUHV— 1 

34 

C 

35 

C 

VARIABLES  BEFORE  VERTICAL  VELOCITY  ARE  MRlTTEN  OUT  THE  SAME 

36 

C 

37 

DO  40  J=1,NUHV1 

38 

- 

IPN (J)=IPR(J» 

34 

00  40  1  =  1  ,NR 

40 

INOl=(I-l)#LR 

41 

XN02=(  1-1ULN 

42 

40 

VW(IN02«Jt  =  VR(lN0l«Jl 

43 

C 

44 

c 

VARIABLE  W=V-WX  IS  NRITTEN  OUT  JUST  BEFORE  ThE  VERTICAL  VELOCITY 

45 

c 

VARIABLE  NX  IS  WRITTEN  OUT  AT  THE  ENO  OF  THE  INPUT  VARIABLES 

46 

c 

47 

00  50  I=I$TART«IN0Ex2 

48 

IN01=(I-1  l*LR  . 

44 

IN02=(T-1I4LW 

50 

VH(IN02«LWI=NZ(II*100. 

41 

50 

VWf  XN02+NUMV  t  =  Vft(  1N01  ♦NUMV  I -Ml  (  It  *100 

c 

VARIABLES  FROM  VERTICAL  VELOCITY  TNR0U6H  THE  LAST  INPUT 

54 

c 

VARIABLE  ARE  WRITTEN  OUT  SHIFTED  OVER  ONE  LOCATION 

55 

c 

56 

NUMV 1=NUMV«1 

«7 

LvlrLW-1 

00  60  J=NUMV 1 tLWl  > 

1PV IJI=IPR(J-1> 

60 

00  60  1=1, NR 

61 

INOt=(  I>1|*LR  •> 

62 

IN 02=1 I-l )*Lw 

*  63 

60 

VN  IIN02*J>=VrI InOI ♦J-j > 

64 

c 

65 

XPlflNUMV  )=*N  • 

66 

JPW( lN )= *WI  • 

i 


45 


87 

C 

88 

C 

SAVE  VALUES  FOR  RETURN  TO  SUBROUTINE  WZcOMP 

89 

c 

- 

70 

WSAVE-WI INR1 

71 

IKD2=(NR-1)*lR 

72 

T  SAVEzVR |IN02+NUMT) 

V3 

AsAy£=VR  IlND2«-NUM  A1 

79 

ULASTS=VRIIN02«NUPV) 

75 

c 

76 

c 

CALL  2 Ml  IT  ONLY  IF  THE  ENTIRE  SE6HENT  HAS  BEEN  COMPLETED 

77 

c 

• 

78 

1  FI  IN0EX2.LT.NRI  RETURN 

79 

CALL  ZMRITIIUNW,1F,01 

80 

RETURN 

81 

END 

Hi  FCHFILE.SIhUl 


OOFMUFCHFlLc'll  .SjHUl  -  . 

C  FUNCTION  SINUL  |N ,A,X,EPS,lNDlc,NRc> 

C  SI 

s  c  when  indic  is  negative*  sinul  Confutes  the  inverse  of  the  m  ? y  si 

9  C  N  MATRIX  A  IV  PLACE.  WHEN  INDIC  IS  2eR0,  SlNUL  COMPUTES  THE  Si 

S  c  N  SOLUTIONS  XI1)...XIn)  corresponding  to  the  set  of  linear  si 

“  6  '  c  EQUATIONS  UITH  AUGMENTED  MATRIX  OF  COEFFICIENTS  In  THE  N  BT  SI 

7  C  NM  ARRAY  A  AND  IN  AOOITION  COMPUTES  THE  I  •  ERSE  OF  The  S} 

8  C  COEFFICIENT  MATRIX  IN  PLACE  AS  ABOVE.  IF  INDIC  IS  POSITIVE ,  SI 

9  C  THE  sEt  OF  LINEAR  EQUATIONS  IS  SOLVED  RUT  THE  INVERSE  TS  NDT  Si 


10  C  COMPUTED  IN  PLACE.  THE  GAUSS-JORD An  COMPLETE  ELIMINATION  HETHODSI 

11  C  IS  EMPLOyEO  WITH  THE  MAXIMUM  PIVOT  STRATEGY.  ROM  AND  COLUlN  SI 

12  C  SUBSCRIPTS  OF  SUCCESSIVE  PIVOT  ELEMENTS  ARE  SAVEO  IN  0»0CR  IN  SI 

13  C  THE  IRON  and  JCOL  ARRAYS  RESPECTI VElV .  K  Is  THE  PIVOT  COUNTER, SI 


19  C  PIVOT  THE  ALGEFRAIC  VALUE  OF  THE  PIVOT  eLe"eNT,  "AX  Si 

is  c  the  number  of  columns  in  a  and  deter  the  determinant  of  the  si 

is  c  coefficient  matrix,  the  solutions  are  computed  in  the  cn*i»th  si 

17  C  COLUMN  OF  A  AND  THEN  UNSCRAMBLED  AND  PUT  IN  PROPER  ORDER  IN  SI 

IB  C  X (1  ). «.X( N|  USING  THE  PIVOT  SUBSCRIPT  INFORMATION  AVAILABLE  Si 

19  C  IN  THE  IROW  AND  JCOL  ARRAYS.  THE  SIGN  OF  THE  DETERMINANT  IS  Si 


20  C  ADJUSTED*  IF  NECESSARY,  BY  DETERMINING  IF  An  ODD  OR  Even  NJN3ERSI 

21  C  OF  PAIRMlSE  INTERCHANGES  IS  REQUIRED  TO  PUT  THE  ELEMENTS  OF  ThESI 

22  C  JoRO  array  in  ASCENDING  SEQUENCE  where  JORO(lROw<I||  =  JCOL III. SI 

23  C  IF  THE  INVERSE  IS  REQUIRED,  IT  IS  UNSCRAMBLED  IN  PLACE  USING  Si 

29  C  VCl)...ylN>  AS  TEMPORARY  STORAGE.  THE  VaLUe  Of  ThE  DETERMINANTS] 

25  C  IS  RETURNED  AS  THE  VALUE  OF  THE  FUNCTION.  SHOULD  THE  POTENTIALS! 


c  pivot  of  Largest  magnitude  be  smaller  in  magnitude  than  eps*  si 

c  the  matrix  is  considered  to  be  SINGULAR  AND  a  true  ZERO  IS  Si 

48  c  returned  as  the  value  of  the  Function.  si 

29  C 

30  C  REFERENCES  CARNAHAN,  LUTHER  AND  WILKES  119891 

31  C  APPLIEO  NUMERICAL  METHODS.  tflLEv,  NEm  YORK. 

32  C 

33  C  SI 

39  FUNCTION  SIMM.  IN  ,A,X,£PS,IN0IC,NRC» 

35  OIMENSION  IROWISO),  JCOLlSOI,  J0R0I501,  Y  ISO  I,  AINRC«NRCI,  XIN»  Si 

38  C  SI 

37  NAX  =  N  ♦  SI 

38  IF  |  INDIC. GE.O  I  !•«:»*!  SI 

39  C  SI 

90  C  . IS  N  LARGER  THAN  SO .  St 

91  xr  |  N.LE.SO  )  GO  TO  5  $1 

92  URITEI  6  ,2001  *  »I 

93  SINUL  £  0.  SI 


46 


M 

RCTURN 

$1 

45 

C 

SI 

46 

C 

.  BC6IN  ELIMINATION  PROCEDURE  ..... 

SI 

47 

5 

deter  =  1. 

Si 

46 

00  18  K  :  1,  N 

SI 

49 

KMl  =8-1 

SI 

SO 

c 

si 

M 

c 

.....  SEARCH  FOR  THE  PIVOT  ELEMENT  . 

SI 

MlVoT  =  0. 

ST 

00  11  1  =  1,  N 

s) 

54 

00  11  J  =  1,  N 

S! 

55 

c 

.  SCAN  IRON  AnO  JCOL  ARRAYS  FOR  INVALID  PIVOT  SUBSCRIPTS  .. 

...SI 

56 

IF  |  8. £0.1  1  GO  TO  9 

S.I 

<7  ‘ 

00  8  ISCAN  =  1,  8M1 

SI 

00  8  JSC AN  r  1.  KM  1  > 

SI 

IF  1  I.Eq.IroWIIsCAnI  )  60  TO  11 

SI 

60 

8 

IF  (  J.eO.JcOLIjScaNI  )  GO  TO  11 

SI 

61 

9 

IF  1  AOS (A(l »JI ) «LE.  ABS(P1V0T|  )  60  TO  11 

si 

62 

PIVOT  =  A 1 1 .  J) 

SI 

63 

IRON (8 |  =  I 

SI 

64 

JCOL (8  1  =  J 

SI 

6S 

11 

CONTINUE 

SI 

66 

c 

SI 

67 

c 

. INSURE  THAT  SELECTED  PIVOT  IS  LARGE®  TH«N  F  ’S  . 

SI 

68 

IF  1  ABS«PI VOT* «GT»EPS  >  00  TO  13 

SI 

69 

STMUL  =0. 

SI 

TO 

RCTURN 

SI 

Tl 

c 

ft 

SI 

T2 

c 

. update  the  determinant  value .  '■  > 

SI 

73 

13 

iROtf K  z  IRON (8)  . 

SI 

7r 

JCOLK  =  JCOL 18 »  *  * 

SI 

75 

DETeR  =  oeTeR*pivot 

SI 

76 

c 

SI 

77 

c 

.  NORMALIZE  pivot  row  elements  . 

SI 

78 

00  19  J  =  1*  MAX 

SI 

79 

14 

A<lROlfK,Jj  =  A  ( IROtfK  .JJ/PIVOT 

si 

80 

c 

SI 

81 

c 

.  CARRY  OUT  ELIMINATION  AND  DEVELOP  INVERSE  . 

SI 

82 

aiXrohk, jcolk)  -  i. /pivot 

si 

*  ** 

Do  18  I  =  1.  N 

Si 

AIJCK  =  All  *  JCOLK  1 

SI 

«. 

IF  (  I «E  0> I®  OWr  )  GO  TO  18 

Si 

86 

A  II, JCOLK)  =  -  AIJCK /PIVOT 

SI 

87 

00  17  J  =  1,  MAX 

SI 

88 

17 

IF  (  J.Ne.JcOLK  I  Al It  J)  =  A( I v Jl  -  AI JCK*AI lRO«Kt Jl 

SI 

89 

1® 

CONTINUE 

SI 

90 

c 

SI 

9l 

c 

.  ORDER  SOLUTION  VALUES  (IF  ANT)  AN0  CREATE  JORO  ArrAY  .... 

.  Si 

92 

DO  20  I  =  If  N 

SI 

93 

IKONI  =  IRON  ID  - 

SI 

98 

JCOLI  =  JCOL  III  * 

SI 

95 

jord  i  iron  I)  =  jcoli 

SI 

96 

20 

IF  1  InO  IC.GE.O  1  XI  .1C  OLD  =  AllROuI.MAxi 

SI 

97 

c 

SI 

98 

c 

.  ADJUST  SIGN  OF  DETERMINANT  ..... 

SI 

99 

IF(N.EQ.l)  60  TO  24 

100 

INTCH  =  0 

SI 

101 

NM1  =  N  -  1 

SI 

102 

00  22  1  =  It  NMl 

SI 

103 

IP1  =  I  ♦  1 

SI 

104 

00  22  J  =  IPlt  N 

SI 

105 

IF  (  JORolJI  .6Eaj0R0II»)  60  TO  22 

si 

106 

JTEMP  =  J0R0IJ1 

SI 

47 


107 

jCROtJ)  =  JORDCI) 

SI 

108 

JCRD(I)  r  JTLMP 

SI 

intch  =  INTCM  «  I 

si 

22 

CONTINUE 

SI 

111 

IF  <  1  NT CH/2*2»N£ .INTcH 

)  deter  =  -  dfter 

si 

112 

C 

SI 

113 

C 

. IF  inoic  is  positive  return  with  results . 

SI 

1*9 

29 

IF  (  INOIC.LE.O  )  60 

TO  26  . 

SI 

SIMUL  z  DETER 

v 

SI 

. 

return 

r 

si 

117 

C 

SI 

118 

C 

. if  inoic  is  negative  or  zero,  unscramble  the  inverse 

SI 

119 

C 

FIRST  BY  ROUS  ... 

SI 

120 

26 

DO  28  J  =  1,  N 

SI 

121 

DO  27  1  =  1,  N 

, 

SI 

122 

I ROW  I  z  IROW  II  ) 

, 

SI 

123 

JCOLl  =  JCOL(I) 

si 

129 

27 

YCJCOLI)  Z  A(lRoVI.J) 

SI 

125 

DO  28  I  =  1,  N 

SI 

126 

28 

A  (I ,  J)  =  y(I) 

Si 

127 

C 

•  •  •  • 

SI 

128 

DC  30  I  z  l,N 

SI 

129 

DO  29  J  -  1,  N 

Si 

130 

I  ROW  J  z  IROU(J) 

SI 

131 

JC0LJ  z  JCOL(J) 

• 

SI 

132 

29 

Y(lROWJ)  z  A  (I  ,  JCOLJ  ) 

si 

133 

DO  30  J  -  1,  N 

* 

SI 

139 

30 

Afl.J)  =  YfJ) 

' 

SI 

135 

c 

Si 

136 

c 

.  RETURN  For  InDIC 

NEGATIVE  OR  2ER0  ..... 

SI 

137 

SIMUL  z  DETER 

SI 

138 

RETURN 

si 

139 

c 

SI 

* 

c 

.  FORMAT  for  OUTPUT 

STATEMENT  . 

SI 

200 

FORMAT!  10HDN  TOO  BIG  ) 

si 

1*2 

c 

SI 

193 

end 

SI 

PR  T,S  FCHF1lE1.mVFIX1,.VFIX1-s,.CN0S/VFIXI-S,.PRESS/VFIx1-S 

iSV'FCHFILElll  1.VFIX1-S  -  -  -  - 

C  THIS  PROGRAM  COMPUTES  PROFILER  CORRECTED  .  WEK.ICAL 
*  c  velocity  bt  using  the  pressure  derivative 

3  C 

A  C  BEFORE  EXECUTING,  THE  INPUT  D«T«  Feb  Fl*-E  MUST  BE  ALIGNED 

sc  To  some  unit  number,  and  an  output  file  assigned,  time  . 

G  C  INTERPOLATION  of  the  INPUT  FEB  FILE  SHOULD  HAVE  BEEN 

7  C  COMPLETED  !HTP»PR06 .ZlNTERP) 

8  C 

9  C  MAPPING  ELEMENT:  FCHFILE1 .MVFIX1-S 

10  C 

11  DIMENSION  wKSOOl, PUT) 

12  COMMON  KSE6,  IEnO  ,K,K2, 1ND,NV  AR,NSE6,NSS£6  ylUNR  y!HN  ,  ;)UHV 

13  COMMON  /RHOR/LR ,NR,NBR ,NMBR » NMFR » NFR ,NIR*nAR ylpR  _u) 

19  COMMON  /ROATA/VRI20000) 

IS  COMMON  /ROOCF/FOOCRI50)/RDOCI/IDOCR(SO)/RDOCA/ADOCR(ZOO) 

10  COMMON  /OIAGS/MSGR,MSGW,NNNR,NNNV,NNlP(NNFyNNI,V’/,*RST.IV$T 

17  OAjA  NNNRyNNNH  ,NN  IP,NNF,  NNl,  NNA/SCOo  ,Sqqo  ,20  ,50*5 Jt. JO/ 

18  C  .-»• 

19  C  SET  INITIAL  OATA  IX=  NUMBER  OF  POINTS  ON  EACH  SI**  OF 

20  C  THE  PRESSURE  DATA  POINT  TO  BE  USED  IN  COMPU  Nv- THE 

2|  C  DERIVATIVE)  DEL Tz  TIME  INTERVAL  bETVEEN  P'.NTS) 


I 


I. 


r 

i: 

i: 

i: 

i 


22 

C 

23 

,  K  -8 

29 

«1=K*1 

25 

K  2-2  *K  *1 

.. 

2' 

CELT =.0625 

2 

R  HOG  1=. 9955 

28 

CC£FF=3./<K*<K*l»*»28K*l>*0ELT» 

29 

HS6R=0 

30 

HS6H=0 

31 

c 

32 

MRlTE(6t2| 

33 

2 

format  (*  input  unit*  output  unit,  number 

SE6HEL 

S.  START  SE6HENT , * 

39 

«./.•  PRESSURE  VARIABLE  NUMBER.  VERTICAL  velocity 

VARIABLE  NUH3ER*) 

35 

READ (5  ,1 IIUNRtIUNW,NSE6tNSSEG(NVAR .NUMV 

36 

1 

FORMAT!) 

37 

KFLA6=Q 

38 

ISTART=H»1 

39 

CALL  7 RE  AO (I UNR. IF.NSSE6 ) 

90 

IENO=NR 

91 

00  30  KSEG=1.NSE6 

92 

C  ' 

93 

c 

FCR  FIRST  SEGMENT.  CALL  SUBROUTINE  TO  SET 

FIRST 

K  DERIVATIVES  TD  ZERO. 

99 

c 

FOR  LAST  SEGMENT.  CALL  SUBROUTINE  TO  SET 

last  k 

DERIVATIVES  TD  ZERO. 

95 

c 

96 

IF(KSE6.£0.1)CALL  ENOSIWl) 

97 

IF(KSEG.NE.1)ISTART>1 

98 

ifikseg.eo.nsegicall  endsivii 

99 

c 

so 

00  10  JrlSTART.IENO 

5’ 

c 

* 

c 

CALL  SUBROUTINE  TO  READ  SPAN  OF  PRESSURE 

VALUES 

INTO  ARRAY  P 

S. 

c 

59 

call  press(p.kflag) 

55 

KF LAG-1 

56 

SUH=0. 

59 

c 

c 

• 

c 

compute  PRESSURE  DERIVATIVE 

5s 

c 

60 

00  20  II=-K,K 

61 

2o 

SUM=SUM*II9P(K1>II) 

62 

OPOT=COCFF*SUm 

63 

Nl (J)=RHOG1«oPoT 

69 

c 

65 

10 

continue 

66 

c 

67 

c 

CALL  SUBROUTINE  TO  NRITE  NEW  FEB  FILE 

68 

c 

69 

IF  .!KSE6.E0.NSE6)G0  TO  25 

70 

NSSE6-NSSEG-1 

71 

CALL  ZREADtlUNR(IF.NSSEG) 

72 

25 

CALL  WRIT EtWI 1 

73 

WRITE  (6*3  1NSSEG 

79 

3 

format! •  completed  segment  «».i6i 

75 

2F!kSEG.E0.NSCGI60  TO  30 

76 

NSSeG=NSSeg  +  1 

77 

CALL  EREAOiIUNR.IF.NSSEG) 

78 

30 

CONTINUE 

79 

END 

£j*'*rcHFILEl«I  I  .ENDS/VFIXI-S 

C  THIS  SUBROUTINE  SETS  THE  FIRST  K  PRESSURE  DERIVATIVES  TO  7ERO  FDR 
2  C  THe  ri*ST  SEGMENT ,  AND  THE  L*ST  R  EqUAL  TO  ZERO  FOR  THE  LAST 


49 


■”’1 


segment 


SUBROUTINE  ENDS'Ull 
DIMENSION  RICH 

COMMON  KSEG,IEN0,K,K2tlN0,NV*R,NSEStNSSEe,IUNR»IUNWtNUHV 
COMMON  /RHOR/LR.NR.NgR , NHp R , NMpR , NfR ,NlP f NAR  ,IPP  C20I 
COMMON  /RDATA/VPC20000I 

CCMHON  /Rd0CF/F00CRI50I/R00CI/I00CR(5qI/rO0CA/A0OCR(2QO> 
COMMON  /DIAG S/MSG R, MSG W,NNNR .NNNW ,NN IP .NNF .NNI tNN A  »IRST, INST 

I  STOP: NR  -K  ♦! 

11=1  STOP 
I II=NR 

SET  UST  Cycle  NUMBER  FOR  last  SEGMENT 
IFIKSEG.  E0»NSE6)IEND=I  STOP-1 


IFIKSE6. EO.t 111=1 
IFIHSFG.EO.1  I11I  =  N 
DO  10  1=11,1X1 
Hlfl)=0. 

RETURN 

END 


E35  *FCHFILEH1 1  .PRESS/ VpIX  l-S 


THIS  SUBROUTINE  READS  A  SPAN  OF  PRESSURE  INTO  ARRAY  P 


SUBROUTINE  pRES$  I  p.KFL  A6  ) 

DIMENSION  P(ll 

COMMON  KSE6,IENO.K,K2,INO,NVAP,NSE6,NSSEG,1UNR,xUNW,NUMV 
COMMON  /RNDR/lR (NR«nBr ,NMB ft , NHFR t NFR ,NIR ,NAR .1PR t  20) 

COMMON  /RO«TA/VR|?iaaO) 

COMMON  /ROOCF/FOOCRI SO J/ROOC I/IOOCRJ 50 > /ROOC A/AOOcR I  2001 
COMMON  /0IA6S/HS&P,nSGV,NNNR,NNNM,NN1P,NNF,NNI,NNA,IrsT,XwST 

IFCKFLAG.EO.llGO  TO  U 

FOR  FIRST  CALL.  READ  IN  2*KM  PRESSURE  VALUES 

DO  10  M=1,M2 

I  ND-(  m-h*lr»nvar 
P  (Mis  VP 1 1  NO) 

continue 

RFTURN 

FOR  CALLS  OTHER  THAN  THE  FIRST.  SHIFT  PRESSURE  VALUES  DORN  BY  OvE 

Nl=K2-l 
DO  20  M=  l.K  1 
P(M)rP(M«i) 

IN0=1ND+LR 
ICHK  =LP*NR 

IF  At  END  OF  SEGMENT.  CALL  SUBROUTINE  TO  READ  IN  NEXT  SEGMENT 
and  reao  In  last  pressure  value  in  span 

IFIXft0.6T«XCHR|CAlL  NExSEG 
PCK2I=VRI1nD) 

RETURN 

end 


•S’A-  - 


f.S  FCHFIL£l*N£XSEG/VFIXl-s,.VRlTE/VFIXl-S,.MVFIxl-S 


E33 ’^FCHFILElIl  )  .NEXSE6 /VFIX1-S 

C  THIS  SUBROUTINE  READS  THE  NEXT  SE6ICNT  «Np  RESETS  ThE  INDEX 

v  C 


3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

19 


SUBROUTINE  NEXSC6 

C  CNN  ON  KSE6tlENDtXtK2tlND,NV*R,NSEG,NSSEG.IUNR.IUNN 
COMMON  / RHDR/LR ,NR,N6R ,NMPR ,NHFR ,NER ,N1R ,N*R ,IPR |20» 

COMMON  /RDATA/VR I  200001 

COMMON  /ROOCF/FdOcRI SO I/RD0CI/I00CR( SOl/RDOCA/ADoCRI 2007 
COMMON  /DIAGS/MSGR,MSGW.NNNR,NNNW»NNlP,NNF,NNltNNA»IRST,IWST 

nsseg=nsseg«i 

CALL  7  RE  AO ( I UNR , I F,N  SS  EG  ) 

INO-NVAR 

return 

ENO 


E33‘*FCHFILE1I1 1 .WRITE/ VFIX 1-S 


C 

this  subroutine  writes  the  new  feb  file  > 

«. 

c 

f* 

(THE  OUTPUT  FILE  CONTAINS  2  MORE  VARIABLES  THAN  THE  INPUT  FILE! 

3 

9 

V 

SUBROUTINE  WRITE  (Wit 

5 

DIMENSION  Willi 

6 

COMMON  KSE6,IEN0,K,K2tlN0(NVARlNSC6(NSSE6»lUNRtIUNW»NUMV 

7 

COMMON  /RHOR/LR,NRtNBR,NMBR, NMFR.NFR .NJR.NAR .IP? 1201 

8 

Common  /roata/vri 20cqo i 

9 

COMMON  /ROOCF/FDOCRI SO 1/rOOC 1/ IPOCRl Sol /ROoC A/ AOOCR 1 2001 

10 

COMMON  /DIAGS/MSGR,MSGW,NNNR tNNNWtNNlPtNNF,NNI ,NN A , T RST, IWST 

11 

COMMON  /VHOR/LW,NW,NBw .NMBw.NMFu.NFW.NIW.NAW.IPW^OI 

12 

COMMON-  /WOATa/VUI 200001 

13 

COMMON  /wDOcF/fd0CW|5OI/Wd0cI/IOOcWIS01/W0OCA/AD0CW(2OO1 

19 

c 

15 

DO  10  1=  l.NNF 

16 

10 

FDOCW II IxFOoCR III 

17 

00  20  1=1, NNI 

18 

20 

IDOCW  |I)  =  ID0cR  III  ’  . 

19 

00  30  I=1,NN A 

20 

30 

ADOCW IIl=AOOCRlIl 

21 

NW=NR 

22 

NBWSNRR 

23 

nmbw=nmbr 

29 

~  ' 

nmfu=nmfr 

25 

NFw=NFR 

2' 

NIU-NIR 

2 

NAU-NAR 

28 

Lw=LR«2 

29 

c 

• 

30 

c 

VARIABLES  BEFORE  VERTICAL  VELOCITY  ARE  WRITTEN  OUT  THE  SAME 

31 

c 

32 

NUHVJ=NUMV-1 

33 

00  90  J=1,NUMV1 

39 

IPW W 1-IPRl J| 

35 

00  90  1=1  ,NR 

36 

INol=< I-1191R 

37 

IND2=II-lULw 

38 

90 

VWIIN02«J1=VrIIn01«JI 

39 

c 

90 

c 

VARIABLE  W=V-Wl  IS  WRITTEN  OUT  JUST  BEFORE  THE  VERTICAL  VELOCITY. 

91 

c 

variable  hi  i$  written  out  at  the  end  of  the  input  variables. 

92 

c 

*3 

00  SO  I=1,NR 

99 

INO 1=11-1 19LR 

51 


I  NO  2=  U -1 1*1.11 

¥  w  C  IN  D2*L  d  I  =W  I II I  *1  00 

VWt  lN02«NUMV  l  =  VR|  INDl  ♦NUNV  l-¥l  «  II  *100 


45 

46 

47 

48 

49 

50 

51 
c 

s 

S4 

ss 

56 

57 
S 
5^ 
60 
6| 
62 

63 

64 

65 


50 
C 

Nl«V  J-NUMV«1 
L«1=U¥-1 
C 

C  VARIABLES  EBON  VERTICAL  VELOCITV  TW0U6H  The  last  input  variable 
c  ARE  WRITTEN  OuT  SHIFTED  OVER  ONE  LOCATION 

00  60  JsNUMV  1,LU I 
Ipy  CJI-IPRC  J-i  t 
0  0  60  I=1  *NR 

IN01-C I“il*LR 
1N02=CI-H*tw 

60  VW|IN02*J)sVR|xN01«J-II 

c 

IPy  1  NUMV )?*N  • 

ipw«lni=*wi  * 

C^LL  2 UR  IT  HUNMf  iFtO) 

return 

ENO 


S A*L  IB  C  |  >  .UMORTH  00 - ,* - 

1  C 

2  - c*****  ***»♦»»**♦*♦**»♦»*****»*»********»**  ********** 


3 

-4- 


5 

6 
7 

a 

9 


c  program 

C  PURPOSE- 


C 

C 

C« 


_X_ 

C 


UNORTHOG 
I-IO- 


INSTPUMENT  VELOCITIES  TP  PILES  WITH  THE 
-ORIGINAL  VC  LOCI  T I ES -XN--AC  OU  STIC— AXE-S  COOPOS  . 


•  *- 


10  _ COHMON/Kf  TA-/0RTH06(  4»  31  .WI3-I -  - 

11  C  COMMONS  For  zpead 

12  _ C0MM0N/6H0R/LR*NRtNHHtN»RR«NMrR»NFR»NIRrNAR.,IPRClS» 

13  COMMON/U'AT  A/VP  I  15,500 

16 _ COMMON/ ROOCF/FOCCR(  *»0  ) - 

15  COMMON/RDOCI/IOOCRI 2C > 

_ 16 - COmiOM/.NPOCA/-AEOCR4-Sai — - 


17  c 

IP  - CCMMON/OIAG$/HSGPtMSGWT*NKRTNNNV'fNNlPTNNFTKNItVNAtIfiSTflWST - 

19  C 

20  - C - INI  T  TAL  I2E— CONT  POL  HEADER - - 

21  TATA  LR  , NR, NER, NIR, NAP/15, 500, 40, 20  ,  50/ 

_ 22 _ CAtA— NWKR-yN  NIP-,N*:F-,NM  A/6QU-rlS,6G . 20-.5O/ - - - 

23  c  initialize  velocity  orthocohalization  matrix 

24  _ caU-ORThOg  Cl.l-I .tORlMOGC2  tl4^0RTHOGtJ^-W^ORTHOSI4*l  1/ - 

25  *  .534, .02465,-. 67693,. b576C/ 

26  - DAT  A— OR-THOG  11,21  ,0R  THOG  12  ,2  1*0RT  H 0613*2-1  *0R  ThOB  14,21  / - - 

27  t  .3505, .66162, .01352,. 63666/ 

26 _ DAJLA — QJW«a6tX»3.1,OKT-HOC  C  2, ZL, OR  THO  6(  3,34-, ORTHOS  4  4,31/ - 

29  t  .2268, .03198,. 70948*. 6C185/ 

30  .  - - 

31  MS6RS1 

32  _ MSGWri - : - 

33  C9494444********************************************************** 

_34 _ X* _ EN tER-XNP-UI- X NO— OUTPUT— PC  B— PILE  I H F OR HATZON - * - 

35  C*****  ********  ******|6********************************************* 

37  C  WRlTCf  6, 2001 

38  200  FORMAT!/*  EfclERt  NUlNl  ,NSEG]  ,NSSE61'-I  - - 

J4  Rf  AO  15  ,3001  sUINl,HSEGl,NSSEGl 

40 _ 300—  F  CRHAX-ll . ■  - 


WrlTE*6,90wl  NUlNl.NSrci  .NSSEGl 

42  4 CO  FCRMAU*  NUI*1  =  *,12, *  NSEGl  =  ’  ,15.*  NSSEGl  =  %I51 

43  7C0  F  CRH  AT  1*  NUIN1  =  *,12,  •  NSEG2  =  *.IS»*  9SSEG2  =  *»I5# 

44  *  _ — •  HSGIX2  =  *  ,121  ...  - - - 

«5  C 

46  _ C _ atAO-5RCCIFICA7IONS^FOR-OUTPUI-  FILE - 

97  C 

All  KRIU»6,SrOI  .  -  - 

49  «tO  FORMAT  4*  ENTERS  nugut*> 


so 

51 

52 
55 
59 

55 

56 

57  - 

58 

59  - 

60 

*1 - 

62 


Rt ADC5.3C01  MJOUT  - - 

WTITt  (6,  9««9  »  NUOUT 

9599  F CRMAT  I*.  NUO  UT 121 _ 

C 

C*  PECIM  1>* AI I.  iCCTlON  0F  PROGRAM  HERE 


- 08-l$e0-t00P— *-*SS«X-r-tM.EM— 

CALL  2READ<NUIN1 ,IF,LOOP> 

- po-aoeu  i---a-.hb - ; - 

c 


63  - . ♦♦*♦♦♦»»»♦♦»♦»♦♦»»  »»»>»»»♦»»♦»»»»»♦♦»•»»»»»»  6  — 

69  C*  CONVERT  VELOCITIES  BACK  TO  ORIGINAL  ACOUSTTC  COOROS  • 

65  - C********************************************************** - 

66  C 

- 67- - 00-7  SOO-K=  1.-3 -  - 

6?  -  UlKl  =  VR(K*9.I> 

69  -  C - IF(  I. LT.3»-6RITF|  6,7001  )KtH|K) - - 

70  70C1  FORMAT!  •  K  =  * ,1  3 »  *  Wf  K  >  =  *,Sl0.6> 

71  -  2500 - COAUINUg - 

72  00  2600  K  r  1,3 

- 73 - V— = ORTHOS  <1  ,K7 - 

74  00  2601  7=  1,3 

75  - V— =— V+W IJ)  *ORTHOB|  J* 1 , Kl - 

76  2601  CONTINUE 

77  - VR«*9.I*~SG-V - 

78  C  IFU.LT.31  WRITEC6, 70021  I,K,V 

_ 79 _ 7002 - COOHAI4  S  Isi  .1  a,  i  Kr»,l2t» - MU1  =*.610.61 - 

80  2600  CONTINUE 

81  -  2000 - CONTINUE - 

82  CALL  2WR IT< ~NUOUT» IF , 0 1 

«3  -  1000 - CONTINUE - - 

89  STOP 

.  —OS _ EMO - 


8PR  T ,S  V .PREFIX  # 

•A*L  18  (  L4.ISERPL0T2 - x, - 

1  C  PROGRAM  TSERPLOT? 

- 2 - C - 

3  C  JIM  VESA  CSC  SEPT.  1981-- 

— « - C - : - - - 

5  C  THIS  PROGRAM  PRODUCES  PLOTS  OF  ANY  FEB  VARIABLE 

6  - C — VS-E-IThER— T-IME— OR-OYClE — NUMBS#* - 

7  C 

.  8 - C - »»»»»♦»♦»»»♦♦»»♦»»»♦■*»»♦»  »»**#»»*»»»»  *9  »9»»»»»66»6» »♦»» *♦»*♦♦♦♦»» 99 

9  C  THIS  PROGRAM  IS  EXECUTED  By  THE  FOLLOWING* 

—  10 - c - 

11  C  .  SXqT  VE GA*LI B. TSERPLOT 2 

12  - C - 

13  C  FOLLOWED  BY  TWO  DATA  CARDS  FOR  EACH  SUBPLOT. 
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DATA-C  ARD— 1— CO 


-10 - c- 

15  c 

*46 - C - I4J,IA8  6I6.IP.VMAX,V4TP.V*4N-,I&« - 

17  C 

19 - C - WHErL* - 

19  c  iu  logical  unit  number  of  input  file. 

.  -20  — -C - lABS-IS— -  POlMER-0£TEPHINlWG-*&SCISS-6-84L4ieS - 

21  c  o  -  cycle  number 

-■  22 _ C _ t.  -  R  EC  A-U-V  E-TlME .  I  IN  SECCNUSi - 

23  C  IP  POSITION  OF  ORDINATE  VARIABLE  IN  FEB 

_ _ £ _  _ DATA  ARK-Ay-  VR . _ _ _ _ 

25  C  Y pAX  *  -  MAXIMUM  EXPECTED  VALUE  OF  ORDINATE  VARIABLE 

26  _ c _ JTSTP-* — s-.ORdiNaT-C— LAliE  LlND-1  nTcRVAl - 

27  C  YMIN*  -  MINIMUM  EXPECTED  VALUE  OF  ORDINATE  VARIABLE 

_ 28 _ e _ t. value  qF.  ordinate  at  origin! - 

29  C  I DEC  -  DECIMATION 

30  C _ 

31  C  *  -  THESE  VARIABLES  ARE  FLOATING  POINT  RATHER  THAN  INTEGER 

32  C - — - - - 

33  C  ***  FOR  A  PLOT  AGAINST  CYCLE  NUMBER  (IABSlS=0>  *** 

_ 36 _ C _ DAJ JL-C  A  8  D.  -2 . CONXAIUSa - 

35  C 

36  C  NUMSEG  ,1  BE  &I  N  fcCYClN - - - 

37  "  C 

38  C  WHERE: - - - - - 

39  C  Nl'MSCG  -  NUMBER  OF  SEGMENTS  TO  BE  PLOTTED 

_ 40 _ C _ If-E  6-U1— s-_SE5  HE  Mi — NUM.  £f  «-£E_E4ESX-lES*ENI--T-0-a£-J*LI>7I-£O - 

41  C  CYCIN  -  NUMBER  OF  CyCLFS  PER  INCH  OF  PLOT 

43  C  NUMSEG  AND  IBEGTN  ARE  IN TEGERS  ,C YC1 N  IS  FLOATING  POINT 

45  C  ***  FOR  A  PLOT  AGAINST  RELATIVE  TIKE  (IABSISR1)  *•* 

- 46 - C - QA1A  CARO  2  CONTAINS: - - - 

47  C 

48  C.  IS  TSEC.INDSEC,  IPTIME - 

49  C 

£3  ^  WHE£Cs  _ —  _ ... _ _ _ _ _ _ _ _  _  _ _ _ 

51  c  ietsec  -  plot  start  time  tin  relative  seconds! 

52  _ C_. IFOSLCLj-— PLOJ  -£N0— EIKf— t  IH_.ReLAT  I V  E-  -SECONDS  1 - 

53  C  IPTIMC  -  LOCATION  OF  RELATIVE  TIME  VARIABLE  IN 

54  C  - THE  FEB  OAIA  ARRAY  VR - 

55  C 

56  C  ALL  THREE  PARAMETERS  ARE  INTEGERS. - 

57  - C- - - - 

58  C  THE  FOLLOWING  DATA  CARD  SHOULD  FOLLOW  THE  SECOND 

59  - C  —  DATA  -CAP D-oF-T HE— L- A 5 7- SU BPfcO-T-* - 

60  C 

-61 - C - 99 , 0,0 ,0  .,  0.  ,0  0  --  — — ..  —  —  — 

62  C 

53  _  -c  -TO  -INSURE— PROPER-TE RMI NA-TI-ON-OF-  -BOTH -T HE-  PR-OBRAN -  - 

66  c  and  the  plot. 

65  C  - - 

66  C  *****  NOTES  ***** 

- 67 - C - - 

68  C  THIS  PROGRAM  IS  SET  UP  TO  USE  THE  36  INCH  PAPER 

69  C  ON  THE  7ET A  PLOTTER ,  As  A -RESULT  THE  MAXIMUM -NUNBER-  - 

77  C  OF  SUBPLOTS  PER  PLOT  IS  5. 

71  -  C - - - - 

72  C  THIS  PRC6RAH  PRODUCES  AN  INTERMEDIATE  PLOT  FILE  OM 

-73 - 1 - UN  IT  -2  5. - IT  Is  NECESSARY  TO  EXECUTE— ONE— OF— TUt-OlSSPL  A - 

76  C  POSTPROCESSORS  TO  OBTAIN  A  PLOT  TAPE. 

76  C  IT  IS  HI6HLY  RECOMMENDED  THAT  THE  INTERMEDIATE  PLOT  FILE 
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77  C  ON  UNIT  25- BE -COPIEC  TO  ANOTHER  FILE  IMMEDIATELY  FOLLOWING 

78  C  EACH  EXECUTION  TO  gUA«D  AGAINST  LOSS. 

79  - C - 

80  C 

$1  - C — ****4***********»**************************************  *»***•••»»*- 

82  IrOI  TE  (6,1) 

83  . 1 _ FORMA 7  L'  -FoR  INSTRUCTIONS -TERMINATE  RUN  AND-ADD’,  - 

88  «  •  VEGA# LI 8. INST /PLOT 2  * ) 

_ 8  5 _ C - 

86  C  COMMONS  FOR  ZREAO 

87  -  - -COMMON/  RHDR/LR  »NctNBR  *NM8R,NNFR.,NFR-»MIR  ,  NAR-,IP8  (  I  SI  — 

88  COMMON/ROCCF/FDOcFt  SOI/RpOCI/IOOcR(5CI/RnOCA/AOOCR( ITCI 

89  _ Common/ roat  a/vr  i  is*  soul - - 

90  COMMON/ 01 AG S/MSG P ,H SG W» NNNR tNNNW • MN lP»NNr,MNl ,NNA , I RST,IwST 

_ 91 _ C - 

92  COMMON/PLOT /I  lit  1  Pi  VMAx(VSTP|YMiNtIOLC»YORlG 

98  "  c  initialize  feb  reao 

95  - ...  PATA  LR  , NR, NFR,NIR,NAR/2C»5C0, 50,53*103/ . .  . 

96  TATA  NHNR,NNlP,Nf!f,NNI,NNA/SC0.2C»50»53*  103/ 

- 9  7 - XORJ  - 

98  C 

99  C  SET  CP  PLOT  PAGE - 

1G0  CALL  COHPRS 

101  _ CALL  BGNPL  (01 - 

102  CALL  SCKPLX 

_ 103 _ CALL-NO  tfilia - 

108  CALL  FLAT80 

105  .  CALL  PAGE  (72.0,18.31  .  .  _ _  ...  .  _ 

106  CALL  TITLE  (C ,0 , C,D ,0 ,C , 1 .0 , 1 . r I 

1D7  _ CALL  GRAF  (  £.0,1  .G,  1.0, G.O,  l.Q, |  .Ql  ...... . . 

108  CALL  ENDC-N  (01 

1C9 . ,.C _ 

110  2fl  CONTINUE 

111  C  READ  INPUT  PARAMETERS  . .  . 

112  TEAD  (5,10001  IU,IA8SIS ,lP,rMAx,7STP,VMlw,IDEC 

113  _ _ IF  (IU.LGt99l  GO  TO  998 -  _  - 

-118 - c - 5 - 

115  IF  (IABSTS.EO.il  60  TO  50 

•  116 - C - 

117  rEAG  (5,10301  NUPSEG, IBE01N  ,CycIN 

—118 - XAXjS-=N  Um  SC-G/4G-T-CIN  /^00»J - : - 

119  IF  (XAXIS.GT.71.)  60  TO  900 

120  - GO-YO-AU - 

121  C 

122  .  50 - CONTINUE - - 

123  REAO  (5  1 100  0)  IS TSEC, 1N0SEC, IPTIME 

.  —128 - € - 

125  60  .  CONTINUE 

126  •  ■  C - - - - 

127  IF  (IABslS.NE.l)  CALL  CyCLNM  (MUMSF  G.ImEGIN  ,CTClN,X  a*IS) 

128  - IF-dABSlSrEOTl-l-CALL-FELSCC- (ISISECrlNDSE&.miMEl - 

129  Y0RI6=Y0PIG*6.75 

- 130 - 6 0-T0-»G . - . ■■■  ■  -  -■  -  -  ■  ■■ 

131  C 

132  -  900 - WRITE  (-6,9011— RA  XI S - 

133  CO  TO  20 

139  - 901 - FORMAT  (•  XAXIS-TOO  LARGE— >AX7SS',F-7.?,/« - 

135  8  *  MAXIMUM  ALLOWABLE  Sl2t  IS  71  INCHES*,/, 

— 1 36 . .  6  — *■  -GUAPL-OT  WTLV-QT-C IWHW-  -  -  -  - 

137  C 

138  - 998 - CONTINUE - 

139  C 
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1*0  - -NRITE-16,  7771 - 

1*1  777  FORK*!!*  CALL  OONEPL* » 

1*2 - CAU—oOJltEl - 

143  1000  FORMAT  () 

1**  C - 

1*5  tNO 


arer.s  v.croiiMVPt-aja - 

GA4LIB  (11  .CyCLNM  /FLO  12 - - - 

1  SUBROUTINE  CYCLNM  < NUMSEG , IEESIN , CyCIN, X A Xl S> 

2  - - C - - - - - 

3  c  this  subroutine  is  called  from1  tserrlot2  to 

— * - C — RRO-OUCC— ROOTS -OR- ANY  ■  FEB — VARIABLE— Bs  CYCLE  HUMBCR - 

5  C 

6  - C — JIM -  VEGA - CSC - SEPT-, — 1*01 - - 

7  C 

8  -  C  *************** ****************.********  *#••*******•  *****#*#• 

9  C 

_ 10 _ C - COM  MOMS -F- OB-? RE-AO _ 

11  COMMON/ RM0R/LR.NP,NBR,NM*«tNMFRtNFBtNIR,MABtIP7 Cl  51 

.  12  - fOMMOK/ RpOCF/FOOCRC  *0 )/RD(>CI/lDoCR( 2C1/R0OCA/A30CRC SCI - 

13  COMMON/ROATA/VRI lSt57C) 

1*  - -  COMM0N/0IAG<,/MS6P,MSGw,NNNG-^tNNU,NNlP,NNF.,NNIVINA,IRST,IWST 

IS  C 

- 16 - COMMON/  PLOT  /Ill  1 1  p,YmAX  ,  Y  s  TBfYlUM  t  ID  EC  t  TOR  IS - 

17  C 

18  - (II MEM-SI  ON  X  ARAY I  1C0G1  ,Y  AR  AY  1 1-GOO  I - 

19  DIMENSION  LA8I7)  ,LBLX  (3> 

-  20  -  -  -C - - - 

21  C 

- 72 - XN0l*=4 - 

23  NCYCLE=1 

.  2* _ iarlSEGIN _ 

25  1FIRST=1 

26  _ iPLTro _ 

27  LSTSCG=NUKSEG«IbEGIN-1 

_ 28 _ XQECSIQ.EC*  I  .0 - 

29  C 

30  -  20 - CONTINUE - 

31  IF  (IFIR5T.E0.il  GO  TO  25 

32  - IF-  iNCYCLE.LT.NRJ  GO  TO-30 - 

33  25  CALL  ZREAO  jIU,lF,Ie» 

- 3* - nCTCLESI - 

35  IF  II FI  RST.4EO  •  1  •  AND  .IF.  NE  .0  )  60  TO  998 

36  -  .  ...  .  .  ir  HTERH.eC.11-G0-T0-1.03 - 

37  IF  I1F.NE.0!  GO  TO  ISO 

38  - IBS  IB  n - 

39  C 

_ --*0 _ -30 _ CQM-LLMLE _ 

41  IF  (7FTRST.  NE.  1 1  60  TO  50 

42  C-  INITIAL J2E  SUCRLOT  ON  FIRST  RAGS - 

43  TNCCOE  (13, 3000, LgL XI 

44  3000  FORMAT  1* CYCLE  NUMBER M| - 

45  CALL  XINTAX 

46  - CAlX.LA6EL-lIFJU4Rl-,MMBBrfOCCRllOl,JOOCR  H>.1,3R.^LAB,L»LYI — 

4T  CALL  FHY50R  1 1,25,709 IC I 

48  CALL  TITLE  <LA? ,  IwO  ,LBLX,  100, LELY, « ,XAXlS,4.S! --  ...  - 

49  XMINS18E6IN»5UC.>499. 

80  XSTPsCYCIk*  2. 0  - - - 

51  yMAX?MUKStG*50«J.nNlR 

82  _ C  ALL — GR  AF— I  xf  I  N*iS  IE*  XtAX,  YJUlUXS  TPUN  AXi - 


53  C 

56  c  frame  subplot 

ss  xlnotmtxaxi s»r.o 

56  -  CALL-SI KIPT  J-1.  CO. -.75-1- 

—57 - CALL— CO  mnP  I-  4—1  .Ot5-«2S» — 

SO  CALL  CONNPT  4XLNCTH *S .251 

59  - : - CALL-COW4M— tXLTiSTHj-xT-S-L 

60  CALL  CONNPT  4-1. CO,-. 751 

- 6-1 - C - 

62  C 


63 - IFMOT** 


66  C 

65  - 5n - CONTINUE - 

66  C  SET  UP  AXIS  ARRAYS 

68  YARAY|INOXI=VR(IP,NCYCLE» 

69  .  C - - - - - 

70  INDX-INOX+I 

71  - NCVCLE=NCYCLt-*«)EC - 

72  C 

—74 - ir.-AIA>aCB.A6  >.EQ.  LSTSESl  I WMU1 - - 

76  IF  4IOOCR41 ).E0. 11  XtEPN=1 

75  - C - 

76  IF  ClNOx.LT  .1001  >  GO  TO  20 

—77 - C - - 


78 

— 79 - 

80 

81 _ 

82 

83 _ 

86 

- 85 - 

86 

87  . 

88 

89 _ 

*0 

— 91 - 

92 

93  ~ 


lNOXriNDX-1 

- CAL  L  CU  ax.fc-  AXAflA  T  .TAB.  AY  « I  NO  X  ,  0  > - 

1N0X=1 

- IPLT- IPLLU - 

IF  IITERn.nE.H  GO  TO  20 

—C _ 

ioo  continue 

— C - TEBMlNATr  SUBPLOT - 

CALL  CURVE  IXARA Y . Y AR AY , INU X . 0 I 

- CALL.  CNORR — (01 - 

CO  TO  I  COO 

— 99* - WRITE.  4  b.999>— LF - 

999  FORfAT  ( *  ERROR  IN  ZREAD,  IF=*,I6» 

— ioun - conxlmic _ 

RETURN  . 

- EN0 - 


8PRT,S  V.RELSFC/P|.0T2 


BA*LIB41  I.ReLSECVFLOT? - . . . 

1  SUBROUTINE  RELSEC  4  XSTSeC  »lNoSEC .'iPTlNCI 

— -  2 - C - - - - - — 

3  C  tHIS  SUBROUTINE  IS  CALLED  TO  PLOT  ANY  FER 

— 6 - c — variable  -against  relative  tumUwuccl. - - : - 

5  C 

--  6 - C - - - - - 

7  C  JIN  VEGA  CSC  SEPT.  1981 

—  8 - C - 

9  C  •***•*••*•***********••*********•♦••*****•••••••••*••••*•**•••••*• 

- 10 - C . . . . . 

11  C  COMMONS  FOR  ZREAD 

12  - C0NM0N/IN4DR/LR»NR»NGR  «NNRR»NWFRtNFR^NIRv  •  '.Ru,Ipt  1ST - 

13  COHNONFROOCF/FOOcRC  50I/R00C1/ID0CR4  50T/RD0C  •  AC  CR4 1001 

- 16 - CONMON/-ROAT  A/VRl  !5»  500T - — - : —  •  - 

15  cONNON/OI AG  S/HSGR»N  SGU»NNNR  tNNNN  tNNlPtNNF tNNl  VRA.tRST.lMST 

46 - c -  — : - ■ _ 

17  C0MM0N/PL0T/IUtZPfYMAXtYSTPtYNINtXbLc«Y0RX6 
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DIMENSION  X A? AVI  1000) ,Y  ARAY 1 1030) 

-OXH{NSI0M-l.ARm  tULX<3> - 

WRITE  (6,111*)  IU,IP, YMAX,TSTP,VMIN ,IOEC 


29 

30 

31 

32  . 

33 

— 3* — 

35 

36  . 

37 
36 
3’ 

_*0 — 

*1 

92 

*3 

** 

*5 

.  *6 _ 

*7 

9* 

*9 

50 

51 

52 

53 
S* 

ss 

S6 

S7 — 
Si 

-S9  - 
60 

-61 - 

*2 

-63  — 

69 

65  — 

66 

-67 - 

66 

69 

70 

71  - 

72 

73 — 

79 

75  - 

76 

77  - 


ITER*»=0 

_ KCVC  Le=C _ _ 

ICO UN 1=0 


IFIRSTrl 

- STRTSC=lSTSLC*l. - 

96  WRIT  El  6,991 

99 - FORM ATt,_NUM EER-OP— SIXOWDS  PER  -INCH?*) _ 

RFA0(S,9i)TlKC 

_SLi _ F-CRNAtXl _ 

XLNGTHr|INoSrc-ISTSEC)/TlNC*1.0 
- IF(XLN6TH.6j.71.JwRlT£t6,97) _ _ 

97  ffRHATI*  X  AXIS  TOO  LON*  (mAX=71):  REENTER* ) 

- I r (XLNgTh.&T .71* 160-1096 _ 


18=  IBM 

- CALL  ZRCAO  IIU,IF,IR) - 

ir  1IF.NE.0)  on  To  99 S 
- IF  (IFIRST.EO.O)  60  TO  315 _ 

initialise  subplot  on  first  pass 


3.700  FORMAT  (A6t  *9*1 

..  CALL  XINTAX  ....  ....  .  _ _  .  .... 

CALL  LALEL  tIPRI IP) ,NHBR,FDOCRI 10), 100CR(11 I, IP, LAB ,LBLY1 

...  CALL  PHYSOR  (1.25.Y09IC1 _ _ _ 

C  tRITF  |6,U13|  rPRZS 

C11U _ FORMAT-1/1.3) _ 

CALL  TITLE  lLAn,U0(L8LX,100tL9LY,6,XLN6TH,*.S) 

CALL  Y6RAF  (SIR T Sc*  TI  Nc  «  YMlN,  YS 1 P  ,Y MAX) _ _ 

ifirst=u 

_  FR1=NRtI0CC _ 

-c - 1 - - - 

25  CONTINUE 

- IF — <-rNO  /.(-T-GT-SC-C  .CC  .  Y<H-IP4Wt4  1 .  AN  OrIS«€CTLE  ,¥  R  < I P  T 1  Me,  NRHI- 

'  *  60  TO  20  •’ 

-« - : _ : _ 

30  CONTINUE 

- NC*CLC-NCVCLE*4 - 


IF— USTSCC.6C. WHIP  TIME, WC+CLEl  1-60-70-3^ 


-  CEL  Ta  TrVR (I PTIMf ,Nc YCLC «IOE C) -VR ( 1»TImE .NCYClE) 

- IF  IOEL TAT, OE.S.D-tO  - T0-20C- - 

IF  IVRIlFTZME.NCYCLCI .6E.INDSEC)  GO  TO  100 

-C- - - - 

C  SET  UP  AXIS  ARRAYS 

YARAY«INDXI=YRIIP,NCYCIC) 

-C - WFtT£-H6,ini)VR»IPTlME,NcYCLC),WRIlF,NtYCLE  — 

C  5  XARAY(InOX) ,YARAY(lNOX) 

-Cl  ill - F«RNAT- (9(F7r2-f3X)-) - - - - 


-  ■  — — -  - - 

CALL  25 LAO  (lU.IF.lB) 

-if  iif.ae.o)  oo-to-ioo — 

>cm£=i-IO£C 


continue 

IFIINDX. Co ,1 0004-60 - TO  2 

-CONTINUE - 

WCYCL£=NCYCL£»I0FC 

-TNOJt=TNOXA| - 

CO  TO  5t 


CONTINUE 

■  IT£  NN-I - ; - 

CONTINUE - 

IF  (INDX.Ne  .1000)  lNnX=lNDX-l 

-VtmC-i  t,IU2MJE8H,  INC* - 

TORhAT  «•  CALLlNfi  CURVE  lTfRN=» ,13 .*  IN0X=*.U) 

CALL  CURVE  JX*R*YtYARAY,INPX,0)  ----- -  - 

ICOL'NT-  1  COUNT  *INCx 

-INOX-O  - 

-IF — LITE  KH  .Nf.4 1  —60— TO— TO - 

XFR  A** £2  XL  NS  TM*2  «  5 


raAHf  subplot 

CALL  STKTPT  t-l.C0,-.T5)  - 

- CALL- CO NNPTI-V. 0,5 .2 ST - 

CALL  COnnpT  I XFr AhE ,5.251 

- - CALL-COkNPI-<XFRAHE  ,-.7S) - 

CALL  COKNPT  l-l.CJ,-. 71 )  ' 


CALL  ENDGR  10) 

*0-10-1 000 - 

*RltE.I6,99  9j  Jtf - 

FORMAT  I •  Error  IN  ZrEAD.  IF=»,Ir) 


CONTINUE 
RETURN— 


APRTfS-  V.LACEL/PLOT 


;A*LlBtl  )  .LABEL/PLOT; 


SUBROUTINE  LABEL  IN AN YAR ,NHFR,DAY  ,1 YR.IP »L**»LHY I 


C  THIS  SUBROUTINE  ENCOOES  THE  0AT.4  LITERALS  USfO  FOR 


C  -WIN -VESA - CM - SEPT-r— 19B1 - 

C 

C 

-C — 4UOPLOT— T-IT-LC - 

,  FNCOOE  IRC, 30C0,LAB )  NAMVAR ,NHBR,OA Y,IYR 


59 


C  OfiOZMTE  TITLE 

IF  (IP.CO.l  .0«.IP.EQ.2.0P.IP.EC.^.OR.IP.ro.£.0<». 

- %  iF.co.9.cB^iP-.cc.t»>-t«coof-4^wt-M«eLy>-  w*mua - - - 

IF  (IF* £0*5 .OR • IP.E Q*6.0P»IP*£Q*7*0R»IP*EQ*11*0R, 

- t--lp,£*.l 2-. tM^P.E^.«^6fr.lP.£WI5»-EKCOOEO,3aW  ,L^L  Y>- NRMVaR 

C 

- RET  UR  M - 

3000  FORMAT  (A6f *  CAST  *,A3,» 


T7TT 


•tFR.3.* 


•Ill 


UNCLASSIFIED 
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